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This review is focused on various properties of quantum phase transitions (QPTs) in the Inter- 
acting Boson Model (IBM) of nuclear structure. The model describes collective modes of motions 
in atomic nuclei at low energies in terms of a finite number N of mutually interacting s and d 
bosons. Closely related approaches are applied in molecular physics. In the N — ► oo limit, the 
ground state is a boson condensate that exhibits shape-phase transitions between spherical (I), 
deformed prolate (II), and deformed oblate (III) forms when the interaction strengths are varied. 
Finite- iV precursors of such behavior are verified by robust variations of nuclear properties (nu- 
clear masses, excitation energies, transition probabilities for low lying levels) across the chart of 

r« nuclides. Simultaneously, the model serves as a theoretical laboratory for studying diverse general 

VO features of QPTs in interacting many-body systems, which differ in many respects from lattice 

models of solid-state physics. We outline the most important fields of the present interest: (a) The 
coexistence of first- and second-order phase transitions supports studies related to the microscopic 
origin of the QPT phenomena, (b) The competing quantum phases are characterized by specific 
dynamical symmetries and novel symmetry related approaches are developed to describe also the 
transitional dynamical domains, (c) In some parameter regions, the QPT-like behavior can be 
ascribed also to individual excited states, which is linked to the thermodynamical and classical 
descriptions of the system, (d) The model and its phase structure can be extended in many direc- 
tions: by separating proton and neutron excitations, considering odd-fermion degrees of freedom or 
different particle-hole configurations, by including other types of bosons, higher order interactions, 
and by imposing external rotation. All these aspects of IBM phase transitions are relevant in the 
interpretation of experimental data and important for a fundamental understanding of the QPT 
phenomenon. 
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1 Introduction 

During the last decades, the classical notion of thermodynamical phase transitions [1] has been extended 
in several directions. One of these extensions is connected with finite systems, counting from thousands 
down to just few tens of particle constituents, see e.g. Refs. [2-5]. The corresponding objects are Bose- 
Einstein condensates, atomic clusters, quantum dots and other mesoscopic systems. As the temperature 
or other external parameters are varied, different "phases", i.e. specific structural configurations of these 
systems, suddenly emerge. Finiteness of the system unavoidably results in smoothening of all relevant 
phase transitional observables, but it makes sense to treat these effects as precursors of true phase 
transitions that would take place in an asymptotic regime. In particular, phase transitions are related 
to the scaling of some essential properties with the increasing size. 

Another extension of standard phase transitions is to systems of interacting quantum objects at 
zero temperature. It turns out that when changing the interaction strength in some systems of this 
sort, one crosses a certain critical point of a nonanalytic change (in the infinite size limit) between 
ordered and disordered "phases" , which are represented by two distinct types of the ground state wave 
function. Since at zero temperature thermal fluctuations disappear, the only internal motion that can 
be responsible for the onset of disorder are quantum fluctuations. This situation is therefore referred 
to as the quantum phase transition (QPT). Quantum phase transitions were quickly recognized very 
relevant in a wide range of systems in condensed matter and many body physics, see e.g. Refs. [6-12]. 

Atomic nuclei are at an intersection of both above-outlined extensions of phase transitions. "Classi- 
cal" phase transitions in nuclei are driven by intensive thermodynamical variables — temperature and/or 
rotational frequency — and take several well-known forms: (a) Transition from Fermi liquid to an ideal 
gas of nucleons is observed in multifragmentation phenomena induced by heavy-ion collisions [13-15]. 
(b) Experimental and theoretical studies (see e.g. [16,17]) support the idea that nuclear matter exhibits 
a phase transition from the paired (superfluid) to an unpaired phase [18,19]. (c) An analogous and 
maybe related effect is observed in the lowest rotational band of many nuclei as a sudden growth of the 
moment of inertia [20]. (d) Transitions between different quadrupole deformed shapes are anticipated 
to appear in a wide range of hot rotating nuclei [21]. 

Manifestations of "quantum" phase transitions in nuclei at zero temperature can be mostly observed 
as abrupt changes of nuclear shapes when crossing certain borders in the plane of neutron versus 
proton number, J\f x Z, see e.g. [9,22-24]. Such transitions can be studied theoretically within some 
phenomenological models of nuclear structure, in which the variation of effective interaction strengths 
defines the relevant domain of application on the chart of nuclides, see e.g. Refs. [9,25-28]. 

This review is aimed at quantum phase transitions in the interacting boson model (IBM). Since its 
discovery by Arima and Iachello in 1975 [29-31], this model has played an important role in modeling 
collective motions of atomic nuclei [32-37] . The success of the IBM is based on a simple and elegant 
algebraic formulation, and on a wide span of relevant dynamical regimes. Closely related algebraic 
approaches are applied in molecular physics [37,38], hadronic physics [39,40], and other areas [41,42]. 

The IBM in its srf-boson "incarnation" exhibits both first- and second-order phase transitions be- 
tween spherical, deformed-prolate and deformed-oblate shapes of the ground state [9,43-45]. Although 
the nonanalytic nature of these transitions can be verified only in the unrealistic limit of the infinite 
boson number, N — > oo (while in application to nuclei, N coincides with the number of valence particle 
or hole pairs, so N < 20 in all realistic cases), some finite-size precursors turn out to be very neat 
already at low boson numbers [23]. The shape transitional predictions of the IBM can be applied to 
low-energy spectroscopic data for numerous isotopic/isotonic chains. 

Apart from describing actual data, the IBM can also be considered as a valuable theoretical tool for 
studying some general features of many-body quantum and mesoscopic systems. More than 30 years 



of the model history gradually disclosed that seemingly simple IBM hamiltonians encode a rich variety 
of complex physical phenomena, such as ground- and excited-state quantum phase transitions, effects 
of coexisting regular and chaotic motions, manifestations of various types of exact and approximate 
symmetries etc. The model offers a deeper insight into the roots of these phenomena. In this sense, it 
can be compared to the Hubbard model of solid-state physics, which was designed as a model describing 
rather specific systems, but in the course of time turned into a general testing ground for application 
of diverse theoretical concepts. 

In attempts to understand the origin and consequences of quantum phase transitions in the IBM, 
considerable effort has been spent and a large number of results obtained — see the impressive (but 
certainly still incomplete) list of references [43-117] below. Although these efforts are probably not yet 
completed, this might be about the right time to make a current summary. 

In this review, the focus is set on both experimental and theoretical applications of the interact- 
ing boson model. With regard to the recent reviews of mostly experimental aspects of shape-phase 
transitions in nuclei [23,24], more emphasis is put to the theoretical side of the problem. We start by 
describing some general features of quantum phase transitions in finite bosonic models and by pointing 
out some differences from the lattice models. A considerable part of the review deals with phase- 
transitional features of the simplest version of the model, the so-called IBM-1. It is presented as a solid 
reference frame for experimental studies, which at the same time serves as a theoretical workshop for 
analyzing the QPT underlying mechanisms. Applications in some of the more advanced versions of the 
IBM are illustrated afterwards. At the end, we will discuss recent generalization of the QPT type of 
behavior to excited states. 



2 Quantum phase transitions in finite bosonic systems 

Effective models of nuclear structure contain free parameters (interaction strengths) whose variations 
naturally induce changes of the ground-state wave function. In a finite system such changes are always 
of the "crossover" type, i.e., smooth in all observables. However, in some situations rather sharp 
transitions are encountered, which may signal true quantum phase transitions in the "thermodynamical 
limit". To decide whether this is the case or not, one needs to identify a parameter, let us denote it N, 
that measures the size of the system. Through this paper, the role of K is played by the total number 
of bosons, N. A necessary requirement for calling a given abrupt structural change of the ground state 
a "quantum phase transition" (more precisely, a finite-K precursor of QPT) is that some of the related 
observables become discontinuous in the limit K — > oo. Let us use this feature as a working definition 
of the quantum phase transition in the type of finite many-body systems that we study here. 

Since in this review the focus is set on sharp transitions of nuclear shapes, the collective degrees 
of freedom are of primary interest. These can be represented by bosonic types of excitations, with 
individual bosons occupying a suitable /mite-dimensional single-particle Hilbert space. The theoretical 
framework for this type of description is provided by the family of interacting boson models (Sec. [31), 
whose quantum phase transitions will be studied in the following sections. It turns out, however, that 
in relation to QPTs all finite bosonic models have some general properties in common. We therefore 
start our discussion by summarizing these properties and showing related examples in simpler IBM-like 
models. 
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2.1 Infinite-size limit of finite bosonic models 

Consider a bosonic many-body hamiltonian with one- and two-body (and possibly higher) interaction 
terms conserving the total number of particles: 

H = E + ^J e k b[b k + ^ v Mmn b\b\b m b n + ... (1) 

k k,l,m,n 

Here b\ and bk, respectively, create and annihilate a boson of the fcth type, while e&, and i/a mn represent 
one- and two-body interaction strengths satisfying the condition 

"klmn "ikmn "klnm ^Ikrim <^nmlk <^nmkl "mnlk "mnkl > 

which follows from the exchange symmetry and hermicity. Eq is just an arbitrary energy shift. Let us 
note that (a part of) the subscript values may also be understood as denoting different components 
of the same boson species (e.g. different angular-momentum projections). Afterwards it will turn 
important that one of the bosons can be treated separately from the others, and for this boson we 
reserve the subscript k = 0. We therefore choose the numbering k = 0,1,2, ...,/, where / will be 
shown to coincide with the number of classical degrees of freedom associated with the system. 

The total number of bosons, N = ^2 k b k b k , is conserved by the hamiltonian (ll). Therefore, we 
can consider realizations of the many-body system with different values of N. However, because of the 
fmiteness of the single-particle Hilbert space (its dimension is /+1), the increasing number of interacting 
particles makes the total energy grow faster than linearly with N. Namely, the two-body part will grow 
proportional to N(N — 1), the three-body part (if any) as N(N — l)(iV — 2) etc. Thus the energy is not 
an extensive quantity. To solve this problem, we have to attenuate individual terms by the respective 
factors, which means that we switch to a modified (extensive) hamiltonian H with interaction strengths 
decreasing with N. Expressing then the energy per particle, TC = j^H (an intensive quantity), one has: 

ft = £ + ]v 2_^ Cfc b k b k + jv(JV-I) Z^ v klmn b k b { b m b n + . . . (2) 

k k,l,m,n 

(we keep the adjustable energy shift, now denoted as £q). Strictly speaking, this hamiltonian represents 
a different system than the original hamiltonian (fTl), although via an appropriate readjusting of the 
interaction strengths in Eq. (pi) we can return to the system (pi) for each individual value of N. 
Now, we can introduce self-adjoint coordinate and momentum operators 



q k = ^ (a k b\ + a k b k ) , p k = j= (j3 k b\ - (3* k b k ) , (3) 



or, equivalently, 



b l = ^f- (PkQk ~ ialPk) , h = ^p- (f3 k q k + ia k p k ) , (4) 

with T = a k Pl + a* k j3 k The commutator of q k and pi reads as: 

[qk,Pi]=ij?8ia. (5) 

We want the commutator to be universal, independent of k, so we take a k and f3 k such that T represents 
a constant. For the sake of simplicity we may choose a k — «£ = (3 k — P% — 77s, hence T = 1. 

The value -^ in the commutator n5« plays the role of the Planck constant h. We therefore conclude 
that in finite bosonic systems, in which the scaling H — > TC from Eqs. (pj) and pi) makes the coordinate- 
momentum representation pi) suitable, the iV — > oo limit represents just the classical limit. This is 



rather important finding since in the QPT systems with an infinite single-particle Hilbert space, where 
the scaling (J2J) is not employed, the N — > oo limit represents the transition to the quantum field theory. 
It is obvious that if Eq. Q is substituted to Eq. (J2J), while replacing N(N — 1) by iV 2 for large N, 
we obtain a coordinate-momentum representation of 7i with no explicit dependence on N. Using the 
commutation rule (ph and neglecting the 0{N~ 1 ) contraction terms, the hamiltonian can be written in 
the following simple form: 

W c i(q, p) = e + 1^2 e k (pI + Qk) + \ ^2 Vklmn (PkPlPrnPn + qkQiQmQn) (6) 

k k,l,m,n 



k.l.m.n 



Note that here we have used the above special choice of coefficients a and (3. In fact, as the contractions 
were skipped, Eq. (pi) represents the classical (i.e. N — > oo) limit of hamiltonian pi). 

The classical hamiltonian (pi) can be equivalently obtained using condensate states [8,9,43,44,118- 

123] 



\N,c) oc 



J2°k b 



N 



10) 



Cfc = ^(qk-ipk) 



(7) 



where |0) is the boson vacuum, and taking 



(N c\H\N c) 

W °'"^ = UfrtN.c) =w ^> (8) 

in the limit N -^ oo. The condensate states O respect the conservation of the total boson number. 
Alternatively, one may use the Glauber coherent states [124] 



\(N),d) ocexp 



E rf ^ 



10) 



(9) 



oo. 



both 



which do not have a fixed N, but allow to vary the average (N) = ^ fc |c4| 2 . For N, (N) 
Glauber and condensate states yield the same results if taking c^ oc (N)~ 1 / 2 dk- 

It is easy to show that quantum fluctuations measured by the dispersion of the energy per particle 
in states ([7]) vanish for asymptotic iV: 



(H 2 ) 



\N,c) 



(H) 2 



\N,c) 



0{N~ 



(10) 



This is a consequence of the attenuation of interaction terms in Eq. d2| with an increasing number of 
bosons. As both condensate and coherent states form overcomplete bases in the bosonic Hilbert space, 
the above result illustrates the fact that the asymptotic number of bosons indeed represents the classical 
limit of the present class of systems. In contrast, for an infinite lattice system with a finite range of 
interactions (not hindered by oc iV^ 1 ) a quantity analogous to (10) behaves as 0(1). 



Efc (A + Qk 



There is an additional condition following from the conservation of the total boson number, namely 
= 2, which is equivalent to the normalization (N, c\N, c) = 1. This constraint can 
be used to eliminate one of the degrees of freedom, for instance that connected with the bo boson. 
Indeed, the absolute value |co| can be calculated from the other values |cjt| with k > through the 
normalization condition, while the phase O of c can be chosen arbitrary, e.g. O — =>• Po = 0. We 
therefore obtain a system with / degrees of freedom restricted by the condition J2 k>0 (p 2 + q 2 ) < 2. 
Furthermore, all coordinates qk with k > can be expressed relative to go using the transformation 



Qk l— ► <?fc = v2(?fc/?o)> which is (for p = 0) equivalent to setting Co = 1 in Eqs. (JtJ) and (|9). The phase 
space of the new coordinates (qi, . . . , g&) and the corresponding momenta (pi, . . . ,pk) is unlimited. This 
procedure, however, implies the appearance of characteristic "form factors" in 7i c \, namely the factors 
[1 + J2 k>0 (pl + 1k)]~ rn m the terms with m-body interactions. 

The above results can be viewed as a consequence of a more sophisticated group theoretical proce- 
dure, based on so-called coset spaces [8, 125]. The procedure starts by the identification of a certain 
subalgebra (called maximum stability subalgebra) of the spectrum generating algebra U(/+l) asso- 
ciated with the bosonic system under study (the spectrum generating algebra is formed of bilinear 
products b k bi). In the present case, the maximum stability subalgebra is taken as U(/)®U(1), where 
the single separated degree of freedom is the one connected with the &o boson. The algebraic coherent 
states [125] associated with the factor algebra U(/+l)/[U(f)<8>U(l)] read as 



\N,z) oc exp 



J2 ( z k b l b o + 4 b lh 



fc>0 



bl 



N 



|0>. (11) 



Using the Baker-Campbell-Hausdorf formula, one can prove that Eq. (11) transforms into the form (J7|) 
with c = cosh |z| and c k = m sinh |z|, where |z| = a/X^>o \ z k\ 2 - 

In the following, we implicitly use the above reduced set of / normalized coordinates and the 
associated momenta, but skip tildes from the notation. We utilize shorthand symbols q = (qi, . . . , qf) 
and p = (pi, . . . ,Pf)- An important quantity for the analysis of the ground-state phase transitions is 
the potential energy. It is obtained from the hamiltonian by setting all momenta to zero, hence V(q) = 
7~^ci|p=o- Although the original hamiltonian (|6]) with /+1 degrees of freedom contains only quadratic and 
quartic terms in its potential V(q), the elimination of go creates also cubic terms. Moreover, the factors 
^1 + ^2 k>0 ql)~ m ma ke the Taylor expansion of V(q) infinite. An explicit formula for the potential will 



be given below in relation to the interacting boson model (Subsec. 3.2). Here, the whole procedure was 



described just to elucidate the general method for obtaining V(q) that will be needed in the following. 

2.2 Nonanalytic evolutions with control parameters 

In the classical limit, iV — > oo, the ground state of hamiltonian (fTl) can be identified with the global 
minimum of the potential energy V(q) = T~i c i\ P =o- Coordinates of the global minimum are q m and 
the ground-state energy reads as E = V(q m ). Under "normal" conditions, V and q m depend on the 
hamiltonian parameters e k and v k \ mn in a smooth, analytic way. However, in some cases — and these 
are of interest in this review — the changes of both q m and Eq are nonanalytic for some values of the 
hamiltonian parameters. Such situations are of course well documented in the literature. There exist 
two standard approaches which make it possible to treat nonanalytic evolutions in a systematic way: 
the catastrophe theory and the Landau theory. 

2.2.1 Catastrophe theory 

The catastrophe theory, initiated by Thorn in the mid 1970's, has been broadly advertised among 
all kinds of scientists and technicians. In short, the theory deals with systems in which "continuous 
causes" can lead to "discontinuous effects". Its application in quantum physics was pioneered by 
Gilmore [10], who showed that the catastrophe theory describes and classifies nonanalytic evolutions of 
the ground state properties of some many-body systems. The application of the catastrophe theory in 
the interacting boson model was first outlined by Feng, Gilmore, and Deans [44] and later elaborated 
by Lopez-Moreno and Castanos [47]. 



The key idea opening the whole field is the concept of structural instability [126] . It can be easily 
explained on the potential of the quartic form V(x) = x 4 , which will turn to be very relevant in the 
forthcoming analyses. But consider at first the harmonic-oscillator potential V(x) = x 2 . Add a small 
perturbation to it, V'(x) = x 2 + ef(x), with e being an infinitesimally small number and f(x) a smooth 
function with all derivatives locally (around x ~ 0) restricted by a common bound. Although the 
perturbation distorts the harmonic behavior close to the potential minimum, it is not difficult to see 
that there is no way how it could change the local topology of the problem, i.e., the number and ordering 
of the minima and maxima (if any) around x ~ 0. The harmonic potential is structurally stable. The 
same holds, less trivially, e.g. for a double- well potential V(x) = x 4 — x 2 , with the boundedness of f(x) 
now being imposed within a broader interval including all local extremes of V(x). 

On the other hand, the pure quartic oscillator, V(x) = x 4 , is structurally unstable, since a small 
perturbation ex 2 with e < does change the topology of the problem: the local minimum at x = 



becomes a maximum and emits two minima located symmetrically at x — ±y||e|. If we consider a 
family of potentials 

V(r ] ;x) = (2r ] ~l)x 2 + (l-ri)x i (12) 

depending on parameter r] 6 [0,1], a structurally unstable pure quartic potential oc x 4 is trapped at 
rj = \ in between structurally stable potentials on both sides r\ < | and rj > |. 

Let us consider more general families of potentials with variables x,y,z, . . . (denote their number 
by n) and adjustable parameters a,b,c, . . . (there are r of them). It turns out that for r < 5 all possible 
forms of such families can be smoothly mapped (the transformation affecting the parameters as well 
as the variables) onto 13 "canonical" forms [126]. These forms classify all possible catastrophes in low 



dimension (r < 5). Example (12) belongs to the equivalence class called the cusp. Since ground-state 
phase transitions between spherical and deformed shapes in nuclei and other many-body systems are 
closely related to the cusp, we will introduce this special type of catastrophe. 
The general cusp potential reads as 

Kus P (o, b; x) = x 4 + ax 2 + bx , (13) 

so it has r = 2 and n — 1. In the parameter plane a x b the potential has a single minimum except in 
the cusp-like region 

a<0, |6|< ^VFaf , (14) 

where two local minima exist and one local maximum in between (a bistable form), see Fig. [I] Both 
minima are degenerate at b = 0, so that if b changes from negative to positive values the minima swap 



and the system described by potential (13) undergoes a phase transition. In the region of parameter b 
demarcated by condition (14) the two coexisting minima indicate a phase coexistence interval typical 
for first-order phase transitions. The limiting parameter values for the bistable form of the potential 
are called spinodal and antispinodal points. 

We assume that the system driven through the phase coexistence region is ideally equilibrated, i.e., 
dwells at the bottom of the lowest potential minimum. Crossing the critical point b = (for a < 0) then 
implies that the rate of change of the system's energy flips from the value characteristic for one minimum 
to the value characteristic for the other minimum. The evolution of the ground-state energy has its 
first derivative discontinuous, which is indeed the defining condition of a first-order phase transition. 

Consider now constant 6 = and a decreasing from positive to negative values. In this case, the 
single potential minimum present at a > bifurcates at a = 0, forming two degenerate branches 
characterizing the situation at a < 0. The system can choose any of these branches as they are fully 
equivalent. If doing so, it turns out that the first derivative of the lowest energy varies in a continuous 




Figure 1: The topology of the cusp catastrophe (left). The surface schematically indicates local equilibrium 



points (stable and unstable) of potential (13) depending on control parameters a and b. The bent part of the 



surface demarcates region ([14]) of bistable potentials. Sample potentials from this region are shown on the 
right-hand side. 



way this time, but the second derivative jumps. This is so-called second-order phase transition. Hence 
the cusp catastrophe accommodates both basic archetypes of phase transitional evolution in the systems 
which are of interest in this review. 



2.2.2 Landau theory 

Landau theory of phase transitions was formulated in the late 1930's [127] as an attempt to develop a 
general method of analysis for various types of phase transitions in condensed matter physics (especially 
in crystals). It relies on two basic conditions, namely on (a) the assumption that the free energy is 
an analytic function of a quantity called order parameter, and on (b) the fact that the expression for 
the free energy must obey the symmetries of the system. Condition (a) is further strengthened by 
expressing the free energy as a Taylor series in the order parameter. It is now known that the Landau 
theory — being essentially the mean-field theory — fails in many systems. However, in the class of models 
we look at here it holds. The reason is the above-explained coincidence of the N — > oo and h~ — > limits 
in models with finite single-particle Hilbert spaces and properly scaled strengths of interactions. 

Hamiltonian (fTl) has a number of external parameters. We will now assume a one-dimensional smooth 
path in the multidimensional parameter space, i.e., consider the hamiltonian parameters depending on 
a single real parameter rj. A phase-transitional evolution of the ground state (if any) shows up as a 
nonanalytic change of the ground-state energy E (rj) = V (rj; q m (r])) at a certain critical point r\ = r\ c . 
Since the dependence of the potential V[j]] q) on r\ is smooth (as the path is smooth), the nonanalytic 
evolution of Eq{t]) is always connected to a nonanalyticity in the trajectory q m (^) of the potential 
minimum. We can write 



discontinuity of ^tEq 



^ discontinuity of ^ fc _i q m -v=> kth order phase transition . 



This results from the following sequence of expressions 

^ (&y) '^0 qm 4q m , £E =(£V) +(|:VV) qm -£q m ... (15) 




Aj u \dr/ v I * V v v /qm d^m, dj? 2 -^U \ drj 2y I ' \d v v v Jim dr) 

/ 1 1 . . . .. ^ \ / q ni 

(where the dot represents the scalar product), which can be continued to an arbitrary order. Thus if q m 
jumps, so does Jr-Eo (first-order transition), the jump of ^-q m implies the same for J^-Eq (second-order 
transition) etc. The transition orders introduced here are consistent with the Ehrenfest classification of 
phase transitions. This classification is not applicable in general (since in real systems some derivatives 
may diverge), but it holds within the Landau theory. 

Now we come to the essence of Landau theory [127,128]. At zero temperature, the equilibrium free 
energy coincides with the ground-state energy and the role of thermodynamical variables is taken by the 
hamiltonian external parameters. The phase of the system can be characterized by a suitably chosen 
order parameter £, which in the present context is a certain function of coordinates, £ = £(q). The free 
energy can be expressed as a function of £, we denote it Vh^r];^), and its equilibrium value, obtained 
by minimization, coincides with the ground-state energy: V^r}]^) = E (tj) = V(7/;q m ). Therefore, 
the above-derived properties of phase transitions of different orders hold also after the replacement 

As will turn out later, the order parameter relevant in our case is given by 



e = ± A E^- ( 16 ) 



fc>0 

It represents a radius in the /-dimensional configuration space (the go coordinate was eliminated in the 



way described in Subsec. 2.1) with the sign ±1. Let us note that the existence of phases characterized 



by £ m > and £ m < is an important ingredient of the general Landau theory. The concrete meaning 



of the sign in Eq. (16) will be discussed in Subsec. 3.2 The above definition implies that for £ m = 



the system's ground state is a condensate of b bosons, see Eq. (|7j), while for £ m ^ the ground state 
is represented by a more complicated mixture of more types of bosons. 

We will be mostly interested in transitions between the £ m = and £ m ^ phases. If there is a 
second-order transition of this kind, in its vicinity the order parameter takes arbitrarily small values. 
Therefore, the free energy can be expressed as a power expansion in £, 

vl(v, o = v (v) + a( V ) e + b{ V ) e + c( V ) e + . . . , (17) 

where the linear term was omitted because of symmetry constraints that apply in our case (see Sub- 
sec. 3.2) as well as in Landau's original context. The condition for the second-order phase transition 



between £ m = and £ m 7^ phases at r\ = i] c reads as 

A( Vc ) = 0, B( Vc ) = 0, (18) 

with C{rj c ) > 0. The £ m = phase is located on the A > side, while the £ m 7^ phase on the A < 
side of r/ c . The evolution of £ m on the A < side close to A = is £ m oc yp4|, hence for a linear 
dependence on i] one has £ m oc y/\rj — r] c \. This implies the critical exponent for the order parameteirl 
equal to \. Note that this particular value is specific for the mean-field description of phase transitions 
of the above type [1] and therefore applies in the finite bosonic models studied here. 



1 General definition [1] of the critical exponent A for a quantity A(e), where e = v Vc , is: A = rim e _>o " ln e e ■ 
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On the other hand, the first-order phase transition between £ m = and £ m 7^ phases takes place 
if conditions 

v L (vc,O = v L (vc;0), 1^(7^0 = (19) 

are fulfilled for a certain value £ = £ m 7^ 0. Close to the second-order phase transition, where the 
terms of Vl with £ 5 and higher can be neglected and C > 0, one finds that a simultaneous solution of 



conditions (19) exists if 

B 2 ( Vc ) = AA( Vc )C( Vc ) . (20) 



If higher than quartic terms in the potential (|17|) are neglected, it can be show that T - b 



AC 



represents the only essential parameter of the potential with A,Bj^0 and C > 0. The remaining 



three "unimportant" parameters can be associated with an energy shift and two scale factors. Eq. ( 20 ) 
can therefore be written simply as % = 4. The £ m = phase is located on the T < T c side and 
the £ m 7^ phase on the T > T c side. The spinodal and antispinodal points, demarcating the phase- 
coexistence region, are given by % = ^ = 3.555... and T a = +00. Let us stress that for an infinite-order 



potential (17) these numerical results may in general be valid only in a vicinity of the second-order 
phase transition (although for some infinite-order potentials, including those relevant here, they hold 
everywhere). 



A closer inspection of Eq. (17) reveals that besides the first- and second-order transitions between 
£ m = and £ m 7^ phases, there exists also a first-order transition between £ m > and £ m < phases. 
The critical point for this transition is given by 

A( Vc )<0, B( Vc ) = 0; (21) 

at this point the equilibrium order parameter apparently changes the sign. The £ m > phase exists on 
the B{rj) < side and the £ m < phase on the B{rf) > side. 



Therefore, the thermodynamical potential (17) describes a system with three phases: £ m = 



(phase I), £ m > (phase II), and £ m < (phase III). We see that the first-order phase transition 



between these phases appear on places determined by a single sharp constraint, either Eq. (20) or 



(21), while the second-order transition is limited to places determined by two constraints, Eq. p8| ). 
The second-order constraints simultaneously fulfill both first-order constraints. Therefore, in a general 
parameter space of dimension p the first-order phase boundaries form two hypersurfaces of dimension 
(p — 1) and the second-order phase transition lies in their (p — 2) dimensional intersection. This phase 
structure will be illustrated by the interacting boson model. 



There is a relation between the Landau potential (17) with terms up to £ 4 and the cusp potential 
(13). Clearly, the potential Vl oc £ 4 at the second-order critical point coincides with the germ of the 



cusp catastrophe. The transition between the two forms can be achieved by a smooth transformation 
which does not affect the cusp topology. If setting C = 1 in the Landau potential (the scale), the 
transformation involves just appropriate parameter-dependent shifts x i— > x — Xq and V >— > V — Vq in 



the cusp potential. The mapping (A, B) 1— ► (a, b) between topologically equivalent forms (17) and (13) 
then reads as [103]: 

a = A-lB\ b=\B(\B 2 -A). (22) 

We may therefore conclude that the phase structure sketched in the previous paragraph belongs to the 
cusp equivalence class. 

2.3 Hamiltonians with a linear dependence on the control parameter 

The classical analysis of the previous subsection needs to be connected with specific quantum signatures. 
To this end, consider a class of quantum models with a linear dependence on a single real control 

11 



parameter 77. The hamiltonian reads as 

H( V ) = H + V V, (23) 

where H and V are mutually incompatible terms, [H , V] 7^ 0. For the Landau analysis the linearity 



means that coefficients in Eq. (17) will be linear functions of 77, which alone would not be an important 
achievement. On the quantum level, however, the linearity represents a considerable advantage. In 
particular, the knowledge of the whole spectrum Ei(i]) (i = 0, 1, ... ,n — 1) at a single value of the 
control parameter, together with the complete set of instantaneous matrix elements (ifii(T])\V\' l l ; j(' r ])) > 
where \^iij])) is the ith. eigenstate of H(r]), determine all spectral observables for any value of the 
control parameter. 

In the many-body case, the linear ansatz is naturally satisfied if the control parameter represents 
a weight factor at a certain interaction term of the hamiltonian. Otherwise it can be justified by the 
Taylor expansion of a general nonlinear hamiltonian around the point of interest. The dimensionless 
parameter 7] may, in principle, vary within the unlimited domain rj G (00, +00), but in practice we 
are always interested in a certain restricted interval 77 G [771,772]- Since we can redefine H such that it 



coincides with H(r]i) and absorb the value (772 — Vi) i n V, the form (23) in its full generality can be 
studied using the constraint 77 G [0, 1]. An equivalent expression is then 

#(r7) = (l-77).ff"(0)+r7#(l) . (24) 

H H +V 

Quantum phase transitions are likely to appear if the limiting hamiltonians H(0) and H(l) correspond 
to two essentially different dynamical modes of the system, e.g. such represented by distinct dynamical 
symmetries (examples given below). 

The evolution of the hamiltonian eigenvectors can be expressed as a unitary transformation in the 
Hilbert space that naturally conserves traces of all operators. This yields very simple predictions for 
bulk properties of the spectrum of hamiltonian (23), namely for the average energy E — - ^27=0 ^i and 



the spread of the spectrum measured by the statistical dispersion (AE) 2 = - Yli=o (^ — E) 2 - While 
the spectrum average behaves linearly with 77, the dispersion is quadratic [98]: 

n~E = TtH + t]TtV, (25) 

n 2 (AE) 2 = [nTiHl - Tr 2 # ] + 2r/ [nTr(H V) - TrH TrV} + rj 2 [nTiV 2 - TrV] . (26) 



The minimum of Eq. ( 26 ) is at 



nTr(H V) - TtH TiV 
Vm ~ nTrV*-Tt*V ' l j 

where the proximity of energy levels induces rapid structural changes of the hamiltonian eigenfunctions 
\ipi(r/)), as can be seen from basic perturbation theory applied to H (77 + 5rj) = H{rf) + 5rjV, e.g. from 
the overlap formula 



[Ek(jj) - EM} 2 
with the squared distance of levels in the denominator. On the other hand, for 77 far away from 77 m the 



i^^-i-W£ ^^ (28) 



second term in Eq. (23) totally prevails and the corresponding wave functions approximately coincide 



with eigenfunctions of V. The spectrum just linearly blows up. 



The most interesting physics happens around the minimum of dispersion (26). This applies also to 



ground-state quantum phase transitions driven by 77. If such a transition exists in the given model, it 
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most likely appears at a critical point t] = r] c that lies somewhere close to r] m . We will assume that our 



choice of Hq and V in Eq. (23) was made so that both ?? m and rj c (if any) are contained in the interval 
rje [0,1]. 

Let us have a closer look on situations when quantum phase transitions related to the structure of 
the ground state can typically take place. First, it is clear that since matrix elements of the hamiltonian 
in an arbitrary fixed basis vary with rj in a smooth (linear) way, any nonanalyticity of the eigenvalue and 
eigenvector evolutions can only occur if the Hilbert space dimension increases asymptotically, n — > oo 
(which for the bosonic models considered here means N — > oo). 

Elementary calculation yields the following expressions for the derivatives of the ground-state energy: 



Since the second derivative in Eq. (29) cannot be positive (E{ > E Vi > 0), the ground-state average 
(V)o = (V'olV'lV'o) never increases. A common situation for the ground-state quantum phase transition 
to appear is when V is a semi-positively definite operator or when some constrains do not allow the 
average (V) to become negative. Then it is likely that the continuously decreasing average (V) 
incidentally reaches zero at a certain point 77 c . If so, r] c represents a critical point of the ground-state 



evolution. Indeed, at this point the first derivative of the ground-state energy in Eq. (29) gets fixed, 



j-Eq = for r] > r] c , and the second derivative jumps to zero, -^E = for rj > r] c . From Eq. ( 2Qh 
we know that everywhere on the right of the critical point there is (ipo\V\^i) (^i\V\^o) = for alEi 
(including i = 0), and this yields {V 2 ) = (ipo\V 2 \ipo) = for 77 > r] c . In other words, the ground state 
becomes an eigenstate of V with zero eigenvalue. 

The situation described above constitutes a second-order ground-state phase transition from a "less 
symmetric", (V) > 0, to a "more symmetric", (V) = (V 2 ) = 0, phase. The notion of symmetry 
is invoked here in relation to the spontaneous symmetry breaking: the "less symmetric" form of the 
ground state usually brakes a certain symmetry that the hamiltonian itself maintains. The ground-state 
average (V)o can be used as a quantum order parameter and related to the classical (mean-field) order 



parameter introduced in Subsec. 2.2.2 



As will be discussed below, in case of the second-order QPT a nonanalytic (for n — > 00) change 
of this parameter is connected with a singular growth of the level density at rj — > r] c and E —>■ E . 
Indeed, as follows from the Pechukas- Yukawa theory [129, 130], the evolution of levels with 77 for a linear 



hamiltonian (23) is analogous to one-dimensional dynamics of a 2D Coulomb gas (however, with the 
product charges also subject to specific variations). To make any of the trajectories nonanalytic, one 
needs to produce an infinite local growth of "charge" at the corresponding place. The Pechukas- Yukawa 



approach and its implications for quantum phase transitions will be discussed in Subsec. |7.2 

For the second-order QPT the order parameter changes continuously, but with a discontinuous 
derivative. In contrast, the first-order QPT involves a discontinuous (for n — > 00) jump of the order 



parameter itself, i.e., the discontinuity of already the first derivative in Eq. (29) at r\ = rj c . In terms 
of the Coulomb gas analogy, this needs locally an infinite "force" , which may be caused by a crossing 
of a pair of levels (or their sharp anticrossing, indistinguishable from a real crossing). Although such 
effects are usually accompanied by an infinite growth of the level density, similar to the second-order 
transition, this is not necessarily the case. Therefore, the mechanisms underlying the first- and second- 
order transitions (a sharp crossing or anticrossing of two levels and a local singularity of the level density, 
respectively) may be interrelated, but in general are different. 
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2.4 Simple example: Lipkin model 

We would like to sketch here QPT properties of the model introduced by Lipkin, Meshkov and Glick in 
1965 [131-133], which is often referred to as the Lipkin model. This will take us very close to the case of 
the interacting boson model that will be opened in the next section. One can find extensive literature 
investigating various properties of the Lipkin model, including its phase transitional behavior, see e.g. 
Refs. [7,10,27,134-138]. 

The model is formulated in terms of pseudospin operators { J z , J+, J_} that form the SU(2) algebra. 
A simplified hamiltonian may be taken for example as follows, 

h>(c) = j z -cM j + + j -) 2 i ( 3 °) 



where ( is a control parameter that should not yet be identified with r\ from Eq. (24), see below, and N 



is related to the spin quantum number j = |iV denoting the chosen SU(2) representation. There exist 



various modifications of Eq. (30) with other quadratic combinations of the pseudospin operators in the 



form V = jj[J 2 + J\ + w(J-J + + J+J-)] containing an additional parameter w. 



The usual interpretation of the SU(2) algebra involved in Eq. (30 ) is in terms of fermionic operators: 



n n n 



Jz — — 2 2_^ a i- a i- + 2 Z-e a i+ ai + ' ^± ~~ Z_^ a i± a iT • (31) 

j=l 8=1 «=1 

Here, a\_ and a^_ create and annihilate fermions on the lower single-particle level, while a i+ and a; + 
do the same for the upper level. Both levels have capacity Q and the total number of fermions, 
N = X^j( a l- a i- + a l+ a i+) — N-. + N + , must satisfy the condition N < Q. The hamiltonian (30) and 
its generalized forms conserve the parity 



n = (-) N + = (-) Jz+j . (32) 



Alternatively, assuming J. = | ^^ a, , the Lipkin hamiltonian can be interpreted as describing interac- 
tions in an infinite array of spin-| particles. As the strength of interactions between fermions on both 
levels — or between individual spins — increases with (, the ground state is tempted to switch from the 
form with (N_,N + ) = (N,0) (all fermions on the lower level, or all spins down) to a "diamagnetic" 
form with both average occupation numbers (N + ) and (N_) nonzero. Indeed, this happens at a certain 
critical point £ c , which will be determined below for a slightly modified hamiltonian by the methods 
outlined in the previous subsections. 

Using the Holstein-Primakoff mapping, 



J z = tfb-j, J + = b ] ^2j-Wb, J_ = y/2j -tfbb, (33) 

one can translate the hamiltonian into the bosonic form. However, this type of bosonic representation 
is not convenient for the present purposes since the total number of bosons b^b is not conserved and 
since the boson interactions are of unlimited order (because of the square root). A simpler alternative 
is to employ the Swinger mapping, 

J z = Wt - s j s) , J + = t ] s , Jl = sh , (34) 



where s\s and t*,t create and annihilate two types of bosons. Since the parity (32) can be expressed 
as n = (—)"*, where n t = t't is the t-boson number operator, it is natural to consider s to be a scalar 
and t a pseudoscalar boson. The original pseudospin algebra of the model is now expressed in terms of 
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the spectrum generating algebra U(2) of the system of s and t bosons with N = s^s + tH = 2j. We 



write hamiltonian (30) in a slightly modified form, 



H{rj) = (1 - rj) [-±{tU + *Ws + sh)] + rjnt = H'(() - (J z + (!-() 



N 
2 ' 



(35) 



where the control parameter ( was replaced by 1— 77 e [0, 1]. The limits rj = and r\ = 1 are characterized 
by 0(2) and U(l) dynamical symmetries, respectively, since the corresponding hamiltonians coincide 
with the Casimir invariants of the 0(2) or U(l) subalgebras of U(2). 

The U(l) case (77 = 1) represents a system of noninteracting bosons, with the s- and t-boson energies 
set to e s = and e* = 1, respectively. The term with (1 — 77) introduces interactions between bosons 
of both types. If these interactions are too weak, the ground state for infinite boson numbers will be 
a pure s-boson condensate. At a certain critical interaction strength, however, the ground state wave 
function flips into a mixed condensate of s and t bosons. The critical strength can be obtained from 



the variational analysis with trial states \N, q) oc (s^ + qt Jf ) N \0), as described in Subsec. 2.1, which leads 
to the following expression for the classical potential energy: 



V(ri;q) 



(5r/ — 4)g 2 + rjq 4 



[1+q 



2^2 



(36) 



Taking into account the Taylor expansion (1 + q 2 ) 2 = 1 — 2q 2 + 3g 4 — 4g 6 + 5g 8 — . . . , we see that the 
potential (36) has the general Landau-like form (17) with the cubic term missing. Consequently, the 



nonanalytic evolution at rj c — I (where the quadratic term in the numerator changes its sign) represents 



a second-order ground state phase transition. For r\ > r? c , the potential in Eq. (36) has a minimum at 
q m = 0, hence the ground state is indeed the pure s condensate, as anticipated. For 77 = (?7 C 
there exist two degenerate minima at 



£j <Vc 



Qi, 



± 



be 



8-5e 



(37) 



±-i/le. The critical 



describing the mixed condensates. For e — > 0, the minima converge to as q m 

exponent for the "order parameter" q m is therefore equal to |. 

The Lipkin hamiltonian can be generalized to get also the first-order phase transitions [44]. The 
only way to do this (while preserving the two-body character of the model) is to sacrifice the parity 



conservation. Indeed, as proposed in Ref. [95], modifying hamiltonian (35) to the form 

H X {Tl) = (1 — 77) [-jjitU + SH + Xt J t)(t j S + SH + X t j t)] + 7777 f , 

which contains parity- violating terms such as tHHs etc., one obtains the potential 

(577 - 4)g 2 - 4x(l ~ y)q 3 + (y + VX 2 - X 2 )<f 
WiV (1 + g 2 ) 2 

Here, the cubic term is already present for x^O and the ground state first-order phase transition takes 



(3c< 



(39) 



place. It turns out that for potential (39) the condition (20) holds exactly regardless of the distance 



from the second-order critical point. The first-order phase transition therefore appears at 

4 + x 2 



Vc 



5 + x 2 



(40) 



In the following section we will see that these results are very close to those obtained within the 
interacting boson model. 
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2.5 Historical and terminological remarks 

As mentioned in Sec.[IJ the term quantum phase transition comes from physics of infinite lattice systems 
of spin-like objects interacting via finite-range interactions. The order-disorder phase transitions in such 
systems at zero temperature are driven by external control parameters and can be related to zero point 
motions, i.e. purely quantum fluctuations, of the lattice constituents. The advent of this kind of physics 
was marked by a pioneering work of Hertz [6] in 1976. At present, the QPT field belongs to one of the 
most rapidly growing branches of condensed matter physics [11, 12]. 

On the other hand, the use of the QPT term in the context of models presently studied [7-10] 
might seem slightly confusing, since — as we saw — the infinite- N limit of such models is just the limit 
of classical physics. However, some of the models that belong to the category "finite" are indeed very 
close to those studied in condensed matter physics. This is most evident for the Lipkin model [131-133], 
which can be cast as a hamiltonian describing an infinite chain of spin- 1 particles interacting by infinite- 
range interactions. A seemingly minor difference from the other lattice models — the infinite range of 
interactions — creates the necessity to damp the interaction constant with increasing iV and therefore 
leads to all consequences following from the convergence to the mean-field description with N — > oo. 
Various forms of the interacting boson model [33], including the one with s and d bosons, is then just 
a natural extension of the Lipkin model from U(2) to higher spectrum generating algebras. As will 
be shown in the forthcoming sections, this extension results in a considerable enrichment of the phase 
structure, allowing also the first-order phase transitions to appear [8,9,43,44]. 

The history of interaction- driven phase transitions, as one may alternatively call such phenomena, 
is however much longer (some historical remarks can be found in Ref. [83]). To our knowledge the 
first author who used the term "phase transition" in this context was Thouless in 1961 [25]. It was in 
connection with what is now known as a "collapse of the random phase approximation" (RPA) in nuclei: 
at some critical value of the hamiltonian control parameter the RPA phonon frequency (determining the 
elementary vibrational mode of the system) drops to zero and becomes imaginary beyond this point. 
This can be considered as a microscopic signature of a sudden structural change of the nuclear ground 
state. 

The Lipkin [131-133] and related pseudospin models created a new wave of interest in similar 
phenomena, see e.g. [27]. The unified language for their description (based on coherent states and the 
catastrophe theory) was developed by Gilmore in the late 1970's [7,8]. Note that Gilmore proposed 
the term "ground-state energy phase transitions", or simply ground-state phase transitions, which we 
sometimes adopt. The field continued growing in the 1980's with the discovery of phase transitions 
in the sd-IBM by Dieperink, Scholten, and Iachello [9,43,44] and also in herefrom inspired so-called 
fermion dynamical symmetry model (FDSM) [139-143]. Since the phases in these models are defined 
by the equilibrium shape of a nucleus in its ground state, the related QPT-like phenomena are often 
called shape-phase transitions. 

The recent increase of interest in this field comes back to the 1990's [22,47,48,144,145]. The 
topic was reopened by several authors pursuing different goals, mainly an analysis of nuclear structure 
evolution [22, 48, 50, 51] and the application of so-called quasi dynamical symmetries [144, 145] and 
critical point symmetries [146-150]. Symmetry (regardless of its concrete incarnation) seems to be a 
unifying theme in a great majority of shape-phase transitional studies in nuclear physics. Whereas in 
the infinite lattice models the quantum phase transition separates ordered and disordered phases of the 
lattice, in nuclear-related models the transition is usually between two specific dynamical symmetries of 
the system, i.e. between two different types of order. In this respect, such transitions can be compared 
to structural phase transitions in solid-state physics. 

In studies of quantum shape-phase transitions in nuclei, the interacting boson model, in its various 
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forms, has attracted the major attention. Apart from a comparison with actual nuclear data [9,46,48], 
the IBM soon became an important testing ground for theoretical investigations of general concepts and 
methods, see e.g. Refs. [47,52,53]. Both these aspects are relevant since the IBM with its rich phase 
structure offers a basic example of a system whose features differ in many respects from the traditional 
QPT systems studied in the context of solid state physics. Therefore, the model may provide essential 
hints for deeper understanding of the QPT physics in general. 

3 The interacting boson approximation 

The Interacting Boson Model (IBM) was proposed in 1975 by Iachello and Arima to describe collective 
excitations of heavy or medium mass atomic nuclei [29-31]. This model combined ingredients of the two 
most successful paradigms used in nuclear physics at that time: the shell model and the geometrical 
(or collective) model [36]. The shell model considers the nucleus as an ensemble of weakly interacting 
fermions occupying single-particle orbits in the nuclear mean field. Despite a considerable truncation of 
the model Hilbert space achieved by activating only the nucleons on valence shells, calculations in heavy 
nuclei away from magic numbers were prohibitively complex. The geometric model attacked the nuclear 
many-body problem from the other side: Heavy nuclei can in some situations be considered as droplets 
of a quantum liquid, with elementary excitations identified with highly correlated collective vibrations 
and rotations. In even-even nuclei, the basic constituents of the model are quadrupole phonons which 
can be represented as bosons with spin and parity l w = 2 + . The collective model has been successful in 
describing certain classes of nuclei away from closed shells. 

The IBM is intermediate between these two complementary approaches in that it connects the 
bosonic behavior of the geometric model to the fermionic nature of the shell model. This is achieved 
using the pairing property of short-range residual interactions. Pairwise coupled nucleons (or holes, 
vacancies left by missing nucleons if the shell is more than half-filled) behave much like bosons. The 
energetically most likely combination of two identical nucleons coupled by a short-range force is the 
one with zero total angular momentum. This can be approximated by an s boson, while the bifermion 
combination with angular momentum 2 maps onto a d boson. 

The original version of the interacting boson model, nowadays abbreviated as IBM-1 [33], is ap- 
plicable to even-even nuclei. The IBM-1 does not separate bosons connected with proton-proton and 
neutron-neutron pairs (this is done in an extended version of the model, the IBM-2 [33] , which is suit- 
able for the description of isovector collective excitations) and does not consider bosons connected with 
mixed proton-neutron pairs (these bosons, increasingly relevant as approaching the M ~ Z nuclei, are 
introduced in more advanced versions of the model, IBM-3 and IBM-4 [32]). The IBM-1 also does 
not consider single-nucleon excitations and their couplings with nucleon pairs (these are treated in the 
interacting boson-fermion model, the IBFM [35], which is applicable in odd nuclei). Some modified 
versions of the IBM also include bosons with other spins and parities, such as g (4 + ), p (1~), and / (3~) 
bosons [33]. Various IBM extensions and their quantum phase transitions will be discussed in Sec. [6| 
while here and in the following two sections we will focus on properties of the IBM-1. 

3.1 Foundations and the algebraic structure 

The interacting boson model, including its simplest version IBM-1, benefits from its transparent alge- 
braic formulation. As indicated above, the model building blocks are s and d bosons that represent 
phenomenological images of + and 2 + pairs of valence nucleons or holes and are also closely related to 
the basic quanta of nuclear collective excitations. Unitary transformations among the six states s'|0) 
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and dJjO), with m = 0,±1,±2, generate the Lie group U(6), which is identified with the spectrum 
generating (dynamical) group of the model. The 36 generators of the associated algebra can be written 
in the form bi m bi> m /, where b 00 = s^ and 62m — ^L- 

Although the separate boson numbers n s and n^ are apparently not conserved by the dynamical 
U(6) algebra, the sum n s + n d = N is. It means that collective states of an even-even nucleus with Np 
valence nucleons (or valence holes) are mapped to the IBM-1 Hilbert space of N = ^Np bosons. This is 
in contrast to the quadrupole phonon model, a bosonized version of the geometric model, where bosons 
directly represent quanta of collective excitations, so that their number varies within one nucleus [151]. 
For instance, in the latter model the ground state of a spherical nucleus is treated as the vacuum, a state 
with N = 0, while in the IBM the same state coincides with a condensate of s bosons with given N > 0. 
In spite of these differences, both models can be connected and rooted in an underlying microscopic 
treatment [123,152-154]. 

An s<i-boson hamiltonian with one- and two-body interactions that conserves the total boson number 
and the total angular momentum has the following general form: 



,(') 



l(o) 



H = E + e d n d + ^ v^, [[b{ x &JJW x [&,, x ^\ q . (41) 

The first term is a constant, which may be included to quantify the nuclear binding energy of the core. 
The second term represents the relevant one-body contributions (the s-boson part can be eliminated 
using the relation n s = N — rid)- The third part corresponds to two-body interaction, the coefficients 
v\ \ v y being related to the interaction reduced matrix elements between normalized two-boson states 
with total angular momentum I. We use the standard notation b\ m = (— )' _m 6/(- m ), it is 5 = s and d m = 
(-l) m cL m , and [A<M x B^}f> = Y Jmim2 {km 1 l 2 m 2 \lm)A { ^B ( t^ where A™ and B^\ respectively, 
are rank li and 1% spherical tensors. 

Taking into account the hermicity of the hamiltonian and its required symmetry under the time 
reversal (reality of coefficients), one finds out that only six real interaction coefficients determine the 
properties of the spectrum (up to the constant shift Eq). The hamiltonian can then be rewritten in the 
following multipole form [33]: 

H = E + e d n d + Cl (L-L) + c 2 (Q"-Q^+c 3 {T^-T^)+c 4 (T^-T^), (42) 

L m = Vl6[S x d]« , (43) 

Q x m = s ] d m + dls + X [d ] x d}V> , (44) 

TL k) = [d j x d]^ , (45) 

where we introduced the scalar product (A® ■ B^) = \f% + i[A® x B®]q . The free parameters now 



read as {E , e^, x, Ci, . . . , c 4 }. Three operators L m (m = 0, ±1) in Eq. (43) define spherical components 
of the angular momentum and generate the physical rotational group 0(3) of the model (the hamiltonian 
is a scalar with respect to this group)jj Similarly, Q^ (m = 0, ±1, ±2) in Eq. (44) represent spherical 



components of the quadrupole operator, determining, e.g. the E2 transition rates. Parameter x, which 
can be chosen within the interval \x\ £ [0, ^], is assumed to have the same value in the hamiltonian 
and in the E2 transition operator (so called consistent-Q formalism [155]). 

Numerical procedures exist to obtain the eigenvalues and eigenvectors of the IBM hamiltonian in 
the general case, but the problem can be solved analytically for some particular choices of parameters. 



2 In this paper, we denote the angular momentum operators and the quantum number associated with the squared 
angular momentum by the same symbol. The rotational algebra and the other orthogonal algebras are referred to as 
O(n), irrespective of whether the determinant is constrained to +1 or not. 



These special cases correspond to dynamical symmetries associated with the algebraic reductions 

f U(5) D 0(5) ) 
U(6) D < SU(3) > D 0(3) . (46) 

{ 0(6) D 0(5) J 

The dynamical symmetries associated with the three chains are named vibrational, U(5), rotational, 
SU(3), and 7-unstable, 0(6). Each of them provides a complete basis for the numerical solution. The 



algebras appearing in Eq. (46) are subalgebras of U(6) generated by operators of the type bj m bi' m i and 
their linear combinations; explicit forms are listed e.g. in Ref. [33]. 

With the subalgebras U(5), 0(5), 0(3), SU(3), and 0(6) there are associated one linear and five 
quadratic Casimir operators. Denoting by C„[G] the nth-order Casimir operator of the group G, the 
general IBM hamiltonian with up to two-body interactions can be written in the following way, 

H = E + aCi[U(5)] + 6C 2 [U(5)] + cC 2 [0(5)] + dC 2 [0(3)] + eC 2 [SU(3)] +/C 2 [0(6)] , (47) 

where we introduced yet another parametrization with {Eq, a, . . . , /}. Explicit expressions for Casimir 
operators as well as for the mapping between the alternative parameter sets are given e.g. in Ref. [33]. 
We point at different conventions used in the literature; below we will implicitly employ the definitions 
listed in Ref. [49]. 



If some of the coefficients in Eq. (47) vanish such that the hamiltonian contains only Casimir 



operators of the subalgebras associated with a single reduction in Eq. (46), the system possesses the 
corresponding dynamical symmetry and the energy eigenvalue problem can be solved analytically. The 
hamiltonian eigenfunctions are determined by the hierarchy characterizing the embedding of irreducible 



representations of groups in chains (46). The eigenfunctions are determined by the corresponding 
irrep labels (these represent conserved quantum numbers) and read as \[N],nd,T,U£.,L) for the U(5) 
dynamical symmetry, | [N], A, /i, Kl, L) for the SU(3), and | [N], a, t, ua, L) for the 0(6). Reduction rules 
specifying which labels are contained in the given irrep can be found in Ref. [33] . Since the eigenvectors 
are independent of the actual values of hamiltonian parameters, they make available clearcut predictions 
for observables like transition rates or transfer amplitudes, which can be used to experimentally verify 
whether a particular dynamical symmetry is present. Analytic expressions for energies are as follows: 

H\[N],n d , t, v A , L) =an d + bn d (n d + 5) + ct(t + 3) + dL(L + 1) U(5) , 

H\[N],X,/j,,K L ,L) = e[A 2 + /i 2 + Au + 3(A + u)] + dL(L + l) SU(3) , (48) 

H\ [N},a, t, z/ A , L) = f a(a + 4) + cr(r + 3) + dL(L + 1) 0(6) . 

Note that the energies do not depend on the so-called missing labels v A and Kl, respectively, that 
characterize a nonunique embedding of 0(3) irreps in 0(5) and SU(3) ones. 

Besides the dynamical symmetries defined above, there exist two additional ones differing from their 



respective counterparts in Eq. (46) by the choice of the relative phase between s and d bosons [156-160]. 



These additional symmetries can be represented by the following chains 



"(^oSts)^ * 3 '- < 49 > 



where SU(3) and 0(6) can be obtained from "standard" SU(3) and 0(6) via the gauge transformation 
s* 1— > e^s* with = 7r and ±|, respectively [other phases are avoided on the basis of the required 
reality of the interaction coefficients in Eq. pi] ), i.e. the time-reversal invariance]. Note that since the 
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U(5) chain a priori separates s and d bosons, it has no counterpart in Eq. (49). Although the additional 
dynamical symmetries are closely related to the original ones (for instance, they yield identical spectra), 
the corresponding hamiltonians written in terms of the original Casimir operators look like if having 
no dynamical symmetry. The ratios of coefficients in expansion (47) needed to fabricate dynamical 



symmetries (49) are as follows: 



b : c : e : / 
a : b : e : / 



2 : 2 : -6 : -1 : 4 
4iV+8 : -4 : : -1 



SU(3) 



(50) 



The existence of additional ("hidden" [157]) dynamical symmetries (49) results from a more general 
property of the IBM-1 hamiltonians, namely the existence of discrete transformations in the model 
parameter space that leave the spectrum unchanged (so called parameter symmetries [159]). 



To conclude the basic overview of the IBM-1 properties, we note that all hamiltonians (47) with 



e = 0, which are transitional between U(5) and 0(6) dynamical symmetries, conserve the 0(5) quantum 
number r, called seniority (a bosonic analog of the fermionic seniority [36]). As a consequence, these 
hamiltonians are still integrable — although the U(5) and 0(6) invariants are no more integrals of motions 
in the transitional regime, their role can be taken by the hamiltonian itself. There is no general 
analytical energy formula valid along the [U(5)-0(6)]dO(5) path but for these hamiltonians the original 
diagonalization problem has been translated to an equivalent problem of solving a certain set of algebraic 
equations (a semi-analytic result) [161-163]. The structure of the spectrum and the form of the wave 
functions exhibit some special features along this transition. In particular, subsets of levels with different 
seniority quantum numbers do not interact with each other and therefore cross with no repulsion as the 
model control parameters vary. 



3.2 Geometrical analysis and the shape-phase structure 

A general method for the determination of the classical limit for hamiltonians of the type (flj was 
described in Subsec. |2.1[ This method can be adapted in the IBM case, yielding a coordinate represen- 



tation of the bosonic many-body hamiltonians (41), and simultaneously a geometrical interpretation of 



the equivalent algebraic hamiltonians (47). The classical limit can be obtained from condensate states 
of the type Q, 



\N,a) 



2\N 



JV!(l + £|a m | 2 ) 



E 



-i N 



™m,U' rn 



10) 



(51) 



(here with the normalization coefficients included), equivalent to the algebraic coherent state (11), or 



from the corresponding Glauber coherent states (J9I). The classical hamiltonian is constructed in analogy 
to Eq. @. 

Coefficients a m in Eq. (51) with m = 0,±1,±2 transform as I = 2 spherical tensors and are, in 



general, complex. To extract coordinates q m and the associated momenta p m , the following prescriptions 
can be used [124]: 



a rn — /jj Wm ' l Vr, 



(In 



1 {(-) m a„ m + a; 



V2 



a, 



Vr, 



^2 il m % Vm\ 1 

f 2 [(-) m a*_ m - a m ] 



(52) 
(53) 



Both q m and p m are again the I = 2 spherical tensors, still complex, but since they satisfy q^ = (—) m q- m 
and p* m = (— ) m p_ m , there are only 5 independent real values for each of them. These define 5 degrees 
of freedom associated with the model. 
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By definition, the classical hamiltonian can only contain scalar combinations of coordinates and 
momenta. There are only two independent scalar combinations of the g's, namely [164] 

[?x?r = ^ 2 , [[gxgp)xg] (0) = -y|/3 3 cos3 7 , (54) 

where we included the parametrization through the well known Bohr quadrupole shape variables (3 
and 7 [165]. As seen, (3 represents a radius in the 5-dimensional configuration space, while 7 can be 
associated with one of the hyperspherical angles. The remaining 3 angles are related to Euler angles 
describing a rotation to the frame where both tensors q m and p m become diagonal (the principal axis 
system). The momenta associated with Euler angles appear in the kinetic energy, but the potential 
energy depends only on (3 and 7. It has the following general form [9,43,44, 119, 120] 



A(3 2 + Bp 3 cos 3 7 + C(3 4 



nPn)= n 1 02^2 ' ( 55 ) 



where the denominator follows from elimination of the degrees of freedom connected with the s boson 



see Sec. 2.1). Coefficients A, B, and C are determined from the general IBM hamiltonian parameters, 



e.g. from a, b, . . . , f in Eq. (47). Such expressions can be found in the literature [33] and we do not list 
them here. 



Potential (55) has the general Landau-like form (17) except it has two variables, (3 and 7. However 



we see that since 7 enters the potential (55) only through the cos 37 dependence in the cubic term, the 



minimization in this variable can be performed separately. The global minimum is either at 7 m = 
(equivalently at ^ or 4^) for B < 0, or at 7 m = | (equivalently at tc or ^) for B > 0. The second 
possibility can be expressed via changing the sign of the corresponding j3 m and simultaneously setting 



7 m = 0. Therefore, the global minimum of the potential (55) can be found using a simplified form with 
7 = and (3 G (—00, +00), which is nothing but the minimization of the Landau potential (17) with 
£ = (3 sign (cos 37). We shall stress, however, that this equivalence is justified only if searching for the 
global minimum of V(/3, 7); any dynamical conclusions based on the above procedure (e.g. identification 
of secondary minima with (3 < etc.) would be strongly misleading since the dependence on 7 is in 
general relevant. 



To separate the "phases" determined by the potential (55), the conditions derived in Subsec. 2.2.2 



can be applied [60]. Let us stress that in the present case the phase separatrices given by Eqs. (20) and 



(21) are valid exactly. The IBM phases can be described as follows [54]: 



I. Phase with (3 m = is the "more symmetric" phase. It appears for A > j|~. This phase is interpreted 

as spherical, since j3 m = implies that the ground state is formed by a condensate of only the s 
bosons, \i/}q) oc (st^O) (up to the leading order terms in N), thus ^(V'ol^dlV'o) = 0. Hamiltonians 
with dynamical symmetry U(5) belong to this phase (we assume e s < e<j) . 

II. Phase with f3 m > 0,7 m = is "less symmetric" and appears for A < JL-, B < 0. It is deformed 

since the ground state is a condensate \ipo) oc (s^ + (3 m d' ) N \0), hence -^(^ol^dl^o) > 0. Promi- 
nent representatives of this phase are hamiltonians with the SU(3) dynamical symmetry. The 
deformation in this phase is conventionally considered prolate. 

III. Phase with (3 m > 0, 7 m = |, that applies for A < JU and B > 0, is also deformed, but interpreted 

as oblate, with \i/j ) oc (s^ — (3 m d ) N \0) and ^(^ol^dl^o) > 0. This phase, which is just a mirror 
conjugate of phase II, can be represented by hamiltonians with dynamical symmetry SU(3). 
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SU(3) 



0(6) 



X<0 



X>0 



Figure 2: (a) A symbolic representation of the parameter space of the interacting boson model with its 
dynamical symmetries. The rear side of the pyramid is equivalent to the front side and any triangle is related 



to all the others by a similarity transformation, (b) The extended Casten triangle of hamiltonian (56) with 



phases I (spherical) , II (prolate) , and III (oblate) . The extended triangle maps the front side of the pyramid in 
panel (a). The I— II and I III first-order phase transitions are connected with narrow coexistence regions while 
the "triple point" at the I-II-III intersection represents a second-order phase transition. 



Note that the dynamical symmetries 0(6) and 0(6) do not represent separate phases, but lie on prolate- 
oblate and oblate-prolate phase separatrices between phases II and III. The IBM-1 parameter space and 
its phase structure are schematically depicted in Fig. [2] 

A brief comment is needed here on the transition between phases II and III. According to the 
Ehrenfest classification, this is the first-order phase transition (the first derivative of the ground-state 
energy jumps) [54]. Nevertheless, it is not of the generic type, since it does not generate the zone of 



coexisting prolate and oblate phases. This difference from the original Landau situation (Subsec. |2.2.2 ) 
is obviously connected with the 7 degree of freedom, whose role in the Landau analysis was explained 
above. Indeed, as approaching the B = line from either side, the secondary "minimum" of the potential 
(55) at 7 = with sign/? = — sign/5 m is only a saddle point. At the critical point, the potential becomes 



totally flat in 7. Another peculiarity is that hamiltonians in parts II and III of the IBM phase diagram 
can be connected by a similarity transformation, which preserves the energy spectrum [159,160]. This 
means that, in a certain sense, the phases II and III are just mirror images of each other. 

3.3 Simplified hamiltonian 

In the previous subsection we saw that the geometrical configuration of the ground state is fully de- 
termined by 2 effective parameters — the ratios of coefficients in Eq. (55). However, the most general 



hamiltonians (42) or (47) have 7 free parameters. It is therefore clear that the phase structure of the 

using a restricted, suitably chosen two-parameter 



IBM-1 can be studied — without the loss of generality 

hamiltonian. In the following, we will work with hamiltonians of the type H 



Qfid - k(Q x ■ Q x ), see 



Ref. [155]. The above form still contains 3 independent parameters (e^, k, and x)> but can be easily 
reduced to a two-parameter form by choosing a specific energy scale (e.g. in units of q). While e^ 
represents a single-particle energy of d bosons (with respect to the single-particle energy of s bosons, 
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which is set to 0), the value of k measures the overall strength of interactions involved in the scalar 
product (Q x ■ Q x ) with given value of x> see Eq. (44). Taking into account the discussions in Sec. [2j we 
rescale the above hamiltonian as follows 



H X (V) = Q--V)M-Q X 'Q*)+Vnd. 



For a fixed x this is the form (24) with 77 e [0, 1] and 



H n 



HQ X -Q" 



N 



V = n d +UQ*-Q 



(56) 



(57) 



Let us note that hamiltonian (56) is analogous to the Lip kin model hamiltonian (38) considered above 



It will turn out that also the phase analysis is essentially the same. 

The case 77 = 1 of hamiltonian (56) clearly corresponds to the U(5) dynamical symmetry, a non- 
interacting ensemble of s and d bosons. An alternative parametrization is therefore the one with 
( = 1 — T]. At ( = 0, the ground state is obviously formed by a condensate of s bosons (the spherical 
phase), but as ( increases the rising interaction strength with d bosons tends to change the nature of 
the ground state into a more complicated mixture of s and d boson states (the deformed phase). This 



mechanism is rather similar to that already discussed for the Lipkin model (Subsec. 2.4). However, the 



present system yields a richer structure since the nature of transitions between the noninteracting and 
interacting phases can be further tuned by changing the quadrupole operator parameter x- 



Comparing the simplified hamiltonian with the one in Eq. (42) we see that only two terms of the 



general form are preserved. However, if rewriting Eq. (56) in the form (47), we get 



H x {rj) 



+ &(1 - v)x (x + f )] Ci[U(5)] + U l ~ V)X (x + f ) C 2 [U(5)] 
iv - 1) (l + jfX + h 2 ) <^[0(5)] + jki 1 ~ V)xix + 2v / 7)C 2 [0(3)] 
+ WV7^ - 1)*G,[SU(3)] + i(l - V ) (l + f 7 x) C 2 [0(6)] 



N ( 



(5J 



(this expansion is valid for definitions of the Casimir invariants given in Ref. [49]), which contains all 
relevant Casimir operators. As we already know, the case rj = 1 (x arbitrary) corresponds to the 



U(5) dynamical symmetry. Expansion (58) can help us to classify the other dynamical symmetries of 

side of the parameter range: the SU(3) dynamical 



hamiltonian (|56|). These are located on the r? 

the 0(6) at (rj, X ) 



symmetry at (77, x) 



(0,-^3 



(0,0), and the SU(3) at (77, x) = (0,+^r) [cf. 



4\ 



Eq @]. 

With the restrictions 77 G [0, 1] and x ^ [ — 2 ' 
be imaged as a rectangle. However, since on the 77 = 1 side x does not appear in the hamiltonian, it is 
more natural to visualize the parameter range as a triangle, see Fig. |2J(b). This is an extended Casten 



the parameter range of hamiltonian (56) can 



triangle [54], which in addition to the standard dynamical symmetries U(5), SU(3), and 0(6) contains 
also the SU(3). The 0(6) symmetry is, however, not present. 

Three comments on the extended triangle are in order: (a) Although the influence of the parameter 
X on the dynamics fades away with 77 — ^ 1 (which justifies the use of triangle instead of a rectangle), 
the E2 transition rates are determined by the quadrupole operator Q x and therefore still depend on 
X even for 77 = 1 [166]. (b) It can be shown that mirror conjugate hamiltonians with x — ~ \x\ an d 
X = +|x| are always connected by a unitary transformation, following from the discrete parameter 
symmetry of the IBM-1 [159,160]. Therefore, the x < and x > halves of the extended Casten 
triangle are dynamically equivalent, although, as shown below, describing rotors with different types of 
deformation (prolate and oblate), (c) We notice in Eq. (58) that on the 77 = side of the triangle the 
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hamiltonians with x^± 2 s ^ m contain a contribution proportional to the sum Ci[U(5)] + C 2 [U(5)]. 
That is somewhat confusing since this case should correspond to the SU(3)-0(6)-SU(3) transitional 
regime. Nevertheless, taking into account the identity 



Ci[U(5)] + C 2 [U(5)] = |C 2 [SU(3)] - 2C 2 [0(6)] + 3C 2 [0(5)] + ^C 2 [SU(3)] - C 2 [0(3)] 



(59) 



one can always eliminate the contribution of U(5) invariants on the 77 = side by including the SU(3) 
Casimir operator into the expansion (158]) 



The geometrical analysis (Subsec. |3.2[) of hamiltonian (56) yields the following results [52]: The 



coefficients in the potential energy formula (55) read as 



A(j]) = 577 — 4 



B x {r)) 



Wf X (l-77), C*( V )=r ] -h 2 (l-v)- 



(60) 



Notice that the latter two coefficients depend on both 77 and \. The spinodal point for the deformed- 
to-spherical evolution appears at 77 = i] s = |. At this point, A changes from negative to positive values 



and the potential (55) develops a local (for x/0) minimum at j3 — 0. The critical point is located at 



Vc 



4+fx 2 
5+fx 2 



(61) 



Here the depths of the /3 = and /3 7^ minima become equal. This condition is fully equivalent 



to Eq. (40) already derived for the analogous Lipkin hamiltonian. The antispinodal point, where the 



(3 7^ minimum disappears, follows shortly after the critical point, 77 a (x) > 77*, the exac t location being 
determined by a cubic equation not given here. For x 7^ 0, the interval from i] s to 77 a demarcates the 
phase-coexistence region for a first-order phase transition between phases I and II (x < 0) or I and III 
(X > 0). ' 

For x = (that is B = 0) all the three points — spinodal, critical, and antispinodal — merge at 



Vc(0) 



and the QPT between spherical and deformed phases is of the second order. The critical 



exponent for the order parameter (3 is |. When crossing the line x = in a transverse direction 
within the interval 77 G [0, |], the potential becomes 7-flat and the sign of B changes. It means that 
the (/?, 7) = (/3 m , 0) global minimum (valid at x < 0) moves to (/?,7) = (/3 m , |) (valid at x > 0)- 
Consequently, the line segment x — 0, 77 G [0, |] represents a transition from prolate (x < 0) to oblate 



(x > 0) shapes (phases II and III, respectively). As discussed in Subsec. 3.2, this is a first-order phase 
transition, but with no separate spinodal and antispinodal points. The point (r],x) — (|;0) can be 
interpreted as a "triple point" of nuclear deformations [60] since three phases — spherical, prolate, and 
oblate — join there, see Fig. [2] Such a simple phase structure is in agreement with the Landau theory 
(see Subsec. 2.2.2). 

Let us finally focus on an alternative parametrization of the 0(6)-U(5) trajectory with x = 0, which 
is often used in the literature, see e.g. [70,71]. It reads as 

H'(r ] ') = r ] 'n d + (l-r ] %^P, 

where P^ = ^((P-d'—s*s') is a boson pair operator. In this case the critical point appears at 77' = 



(62) 



hamiltonian (62) is closely related to the form (56) with x = as one can write 



1 

2' 



H°(r)) = V n d + (l- 77) Aptp +(1 _ 77)£C 2 [0(5)] + (1 - 77) (JV + 4) . 



(63) 



(4-377)H' 



4-3-q 
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Figure 3: The dimension of the subspace of IBM states with L ■ 
dots correspond to actual dimensions calculated for N = 10,20, 
formula n « ^iV(iV + 1). 



as a function of the boson number. The 
. . , 300, the line represents an approximate 



Since the Casimir invariant C2[0(5)] = n^fid + 3) — (d* ■ S)(d ■ d) is conserved along the \ = path, 
it takes values t(t + 3), where r is the seniority quantum number. This means that the hamiltonians 
H' Q (r}') and H°{rj) yield identical eigenfunctions and spectra connected just by simple transformations 
following from Eq. (J63l . 



3.4 Computational aspects 

In numerical investigations of the IBM phase transitions, it is essential to calculate various spectroscopic 
signatures up to very large boson numbers. This may present a computational problem, since the 
dimension of the relevant Hilbert space increases. An explicit formula for the total number of states 
with all angular momenta L 



0,2,3,4, 



, 2N reads as follows, 



"-tot 



(AT + 5)! 
5!iV! 



N 5 
120 ' 



(64) 



where the approximation on the right-hand side holds for very large N. Dimensions of individual L 
subspaces can also be determined. This can be done exactly by employing the well known reduction rules 
for the embeddings of various irreps corresponding to algebras in any of the IBM dynamical symmetry 
chains, all ending at the 0(3) invariant symmetry subalgebra. The result for L = is shown in Fig. |3l 
together with an explicit asymptotic formula n ~ j^N 2 . In order to determine the ground-state energy 
and wave function for a finite N, a numerical diagonalization of the IBM hamiltonian has to be carried 
out within this subspace. We see that although its dimension grows only quadratically with N, the 
boson numbers N ~ 100 may already cause some troubles. 

In fact, the most difficult part of the computation is not the diagonalization procedure itself. Very 
efficient diagonalization algorithms are commonly available, especially for sparse matrices (thus applica- 
ble in the IBM case), which make it possible to go up to very high dimensions. The part that presently 
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seems to set the limits is the initial calculation of the hamiltonian matrix to be diagonalized. In a 
majority of codes, the evaluation of matrix elements relies on the U(5) dynamical symmetry basis and 
makes use of the coefficients of fractional parentage [36]. These are computed in a recurrent way, the 
complexity growing very quickly with increasing dimension. An alternative scheme using the 0(6) basis 
has been developed [172], but it does not solve the problem of large boson numbers. 

Numerical difficulties can be overcome in the case of the transition between U(5) and 0(6) dynamical 
symmetries along the path conserving the underlying 0(5) dynamical symmetry. As mentioned above, 
transitional hamiltonians of this type are fully integrable. The Casimir invariant of the 0(5) subalgebra 
represents an additional integral of motions, yielding the conserved seniority quantum number r. The 
dimension of the L = r = subspace of states (including the ground state) grows asn~ ^N, which 
is a much slower increase than that in Fig. [3j An additional advantage is that the hamiltonian matrix 
in this subspace can be calculated directly, using explicit expressions for matrix elements of the 0(6) 
Casimir invariant in the U(5) basis given by Arima and Iachello [31]. This makes it possible to increase 
the upper bound for boson numbers by several orders of magnitude (cf. Ref. [108]). 

Some alternative, more sophisticated methods of solutions of the eigenvalue problem along the [U(5)- 
0(6)]dO(5) transition have been developed. One of these methods is based on an infinite dimensional 
algebraic technique. In the IBM context, the method was introduced by Pan and Draayer [161] and by 
Dukelsky and Pittel [162]. The general solution resorts from an exact treatment of the pairing model 
invented by Richardson [167, 168] in the 1960's — see Ref. [163] for an overview and more references. 



If the U(5)-0(6) transitional IBM hamiltonian is written in the pairing form (62), the diagonalization 
problem for seniority r is transformed to the problem of solving a coupled set of n = | (N — r) nonlinear 
algebraic equations (Richardson equations). For large boson numbers N and small values of r it is still 
difficult to obtain all solutions, but to determine the lowest eigenvalues is feasible. The method was 
applied in a QPT analysis by Arias, Dukelsky, and Garcia- Ramos [66]. 

Another approach to the [U(5)-0(6)]dO(5) transition makes use of the so called continuous unitary 
transformation technique, introduced in the 1990's by Glazek and Wilson [169, 170] and independently 
by Wegner [171] — for a compact overview see Ref. [90]. The unitary transformation U that diagonalizes 
the given hamiltonian H is found as a limit U = liniA^oo £^(A), where the dependence U(\) can be 
viewed as a continuous flow in the operator space from an arbitrary initial point U(0). Equations 
governing this flow in such a way that it converges to the solution of the eigenvalue problem have been 
derived and they can be iterated numerically. In the IBM, the method was for the first time employed 
by Dusuel, Vidal, Arias, Dukelsky, and Garcia- Ramos [90]. 

A promising tool for calculations with very large boson numbers was recently presented by Ho, 
Rosensteel, and Rowe [105]. It is based on a generalized equation-of- motion approach put into the 
algebraic framework; for an explanation and original references see Refs. [105,173]. The method com- 
pletely avoids the diagonalization step of the computation and proceeds directly to the construction of 
a matrix representation of the dynamical algebra of the problem. This makes it possible to efficiently 
approximate the relevant observables — energies and matrix elements involving low-lying states — using 
matrices whose dimensions are vastly reduced compared to the full Hilbert space dimension. This tech- 
nique has been so far applied in the second-order (7-soft) shape phase transition within the geometric 
model [105] and in the Lipkin model [173], both cases being analogous to the U(5)-0(6) transition in 
the IBM. 
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4 Quantum phase transitions in IBM-1: theoretical results 

Basic results of the mean-field analysis of ground-state shape-phase transitions in the interacting boson 
model-1 were already discussed in the preceding sections. The use of coherent states for the study 
of IBM-1 ground-state transitions (a technique introduced by Gilmore [7,8]) was for the first time 
reported by Dieperink, Scholten, and Iachello in 1980 [9]. Interpretation of these results in terms of 
the catastrophe theory was given by Feng, Gilmore, and Deans in 1981 [44] and later, in 1996, by 
Lopez-Moreno and Castaiios in a detailed analysis [47]. 

The latter work also noticed the phase transition between prolate and oblate shapes (the transition 
achieved by changing the relative phase between s and d bosons). This was then independently discussed 
in Ref. [54], where the extended Casten triangle of Fig. [21(b) was introduced. The IBM-1 shape-phase 
structure with shape types I (spherical), II (prolate) and III (oblate) was put into the context of Landau 
theory by Jolie, Cejnar et al. in 2002 [60,67]. This kind of phase diagram, with a "triple point" in the 
middle, represents the simplest possibility allowed by the Landau theory with a single order parameter 



j3 (for the discussion of differences from the original Landau context see Subsec. 3.2). Atomic nuclei 
seem to belong to its few realizations in nature [174]. 

From the theoretical point of view, an attractive property of the IBM is the simultaneous presence of 
both first- and second-order phase transitions. This naturally allows for comparative analyses probing 
the general features of QPTs. As foreseen in Ref. [44] and newly proposed by Vidal et al. [95], the 
same structure can be reproduced also with an extended Lipkin hamiltonian containing parity violating 



terms (Subsec. 2.4). However, if one sticks to models with two types of bosons — one scalar and one 
non-scalar — under the strict assumption of angular momentum or parity conservation, the IBM with s 
and d bosons is the simplest case with QPTs of both kinds [103]. 

4.1 Calculations beyond the mean field 

So far, the phase transitional properties of the interacting boson model were discussed within the 
mean field (or Hartree-Bose) approximation, i.e. the approach based on the coherent state formalism 



sketched in Subsec. |2.1| This approach yields correct results in the infinite-size limit, N — ► oo, but to 
predict finite-size corrections, one needs to employ more sophisticated techniques. From the theoretical 
viewpoint, this is an extremely important task since the scaling of spectroscopic observables with N 
plays the key role in describing the precursors of quantum phase transitions in finite systems. Indeed, 
in the quantum critical points the scaling becomes singular (different from other points of the model 
parameter space), which may be considered as a defining feature of QPTs in finite systems. 

To calculate the finite-iV corrections, diverse strategies are followed. A direct possibility is to 
perform numerical diagonalization of the IBM hamiltonian for very large boson numbers and to analyze 
the convergence of the relevant quantities of interest to their asymptotic values. In practice, this route is 
viable only for the transition between U(5) and 0(6) dynamical symmetries that conserve (all the way) 
the underlying 0(5) dynamical symmetry, i.e. for the second-order phase transitional path. Indeed, 



as explained in Subsec. |3.1 , the model is integrable along this path, which makes it possible to obtain 



results up to very large boson numbers using various more or less sophisticated techniques (Subsec. 3.4), 
ranging from the Richardson equations [66] to just a "brute force" diagonalization [103]. Depending on 
the actual set of states included in such U(5)-0(6) calculations (i.e., in relation to the specific analysis), 
the upper values of the boson number can be N = 10 4 [85] or even N = 10 5 [103]. On the other 
hand, purely numerical calculations along the IBM first-order phase transitional paths can be currently 
performed with boson numbers not exceeding iV ~ 10 2 [78,83]. 

Another route to large size calculations leads through expansion techniques for various observables 
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Figure 4: The gap A as a function of r\ around the first-order (left) and second-order (right) phase transition 
for the Lipkin hamiltonian (38). The spherical phase is at rj > rj c , with f]c = j\ for x = ~tk an d f] c = | for 
X = 0. Lines and the full dot are analytic results, open dots correspond to a numerical diagonalization. For 
X = both first and second excited states are included. The inset shows the scaling of A at the critical point 
with N for various values of x- Adapted from Vidal et al. [95] . 



in powers of N. The advantage of such methods lies in their enhanced potential to determine the 
asymptotic scaling laws for both first- and second-order phase transitions. An interesting approach 
[89,90,95, 104] is based on the Holstein-Primakoff type of mapping of s and d bosons onto boson b, 



u, v 



btK, 



s's 



N -n b 



4 s 



( s t d )t ,_ bWN-n b 



(65) 



where rib = E«^^- The following step is a shift transformation c^ = bj^ — vAA* with a complex 
vector A M satisfying J^ (AJ 2 = ||A|| 2 G [0, 1]. This implies (n c ) = (n b ) — ||A|| 2 A^ for states with a fixed 
total boson number, and therefore allows for macroscopic average occupation numbers (n&) while at the 
same time (n c ) ^C N. The resulting hamiltonian is written in terms of 5 pairs of operators cL c M and 
the number N, which enables one to naturally separate terms of different orders in N. 

The mean-field approximation represents the highest order term oc A^ 1 . Note that here we consider 



an extensive IBM hamiltonian, such as the one in Eq. (56 ). For hamiltonians in the intensive, i.e. energy 



per-boson form, Eq. ^, the orders must be of course reduced by 1. To proceed beyond the mean field, 
one needs to diagonalize — order by order — the respective higher hamiltonian terms. Calculations of 
this type have been performed by Arias, Dukelsky, Dusuel, Garcia-Ramos, and Vidal for the second- 
order transition [89,90,95,104] as well as for the first-order transition [95,104]. The first nonvanishing 
quantum correction is of order A^° and can be treated with the aid of the Bogoliubov transformation. 
It is relatively straightforward in the spherical phase, while in the deformed phase one needs to proceed 
separately for angular momentum projections \x = 0, ±1, and ±2 (corresponding to the /3, Goldstone, 
and 7 excitations, respectively) [104]. Corrections to the resulting expressions are of order 0(N~ l ). 
Explicit calculations of terms oc N~ k with k = 1, 2 have been performed within the continuous unitary 



transformation formalism (briefly mentioned in Subsec. 3.4) for the second-order transition [89,90]. 

Without going to technical details, we will report here results of these calculations for the gap A 
between the ground state and the first excited state. The convergence A(r/ C ) — > with increasing A^ at 



2s 



the quantum critical point r] c represents an essential QPT precursor. Note that for an IBM extensive 
hamiltonian far away from r] c the average gap is expected to be independent of N, but in the critical 
region the low levels systematically reduce their spacing. Calculations [89, 90, 95, 104] are consistent 
with the scaling laws 

A(r/ C ) oc N~s (second order transition) , A(r] c ) oc e~ a (first order transition) , (66) 

where the value of a depends on the place where we cross the first-order transition, e.g. on the value 



of x i n the parametrization (56). A simple reasoning for the scaling law in the second-order transition 
was given by Rowe [70]. It is based on the fact that in the second-order critical point the mean field 
potential is the pure quartic oscillator V oc Nj3 4 (factor N, which was not present in the notation of 
Sec. 2 is due to conversion to the extensive form). At the same time, kinetic terms T oc —Nh 2 -^ are 



of order N 1 (since h oc N 1 , see Subsec. 2.1). Therefore, after the transformation /3 i — ^ /5 oc iV 1//3 /3 the 
hamiltonian produces an overall scaling factor iV -1 ' 3 , in agreement with Eq. (66). 



The behavior of the gap in both first- and second-order QPTs of the Lipkin model is illustrated 



in Fig. El adapted from Ref. [95]. It was obtained for hamiltonian (38) which yields an equivalent 
phase diagram as the sd-lBM and the same scaling properties, but at the same time enables one to 
perform numerical calculations up to higher boson numbers (N = 5000 in our case). In the second-order 
transition (x = 0), the ground state and first excited state form a nearly degenerate parity doublet for 
7] < rjc, therefore the figure shows two lowest excited states. Scaling properties of various observables 
in a wider class of two-level boson models, including the sd-IBM, can be found in Refs. [89,90, 104]. 

4.2 Thermodynamical analogies 

Structural changes of the ground state in the quantum critical region were investigated by Cejnar, 
Jolie, Zelevinsky and Sokolov [49, 52, 53] by adopting concepts from statistical physics. These analyses 
define suitable entropic quantities, the wave-function and von Neumann entropies, that make it pos- 
sible to follow a parallel with thermodynamics. The wave-function entropy [158] is just the Shannon 
(information) entropy that measures the spread of an arbitrary wave function \ip) = ^2 ~ on \j B ) in basis 
B = {\j B )}j=i t 2,... of the Hilbert space. For the zth eigenstate of the hamiltonian H{rj) this yields: 

wfw=-Ei«Swr in i a S^)i 2 - ( 67 ) 

3 

Since a phase transition, being essentially a sudden restructuralization of the ground state and near 
excited states, can be seen as a transition between two types of basis vectors (those relevant at each side 
of the transition), the wave-function entropy is well suited to demonstrate basic QPT effects [52,67]. 

On the other hand, von Neumann entropy is constructed to be independent of the basis. It can 
be obtained by adding a small noisy component Srj to the model control parameter t], which results in 
turning all the hamiltonian eigenstates \i^i(r])) to density operators Pi(rj), expressing statistical character 
of the place in the parameter space where the eigenstates are evaluated [175] . The zth state von Neumann 
entropy then reads as 

S i (r t ) = r rr[-p i (r i )\np i (ri)], (68) 

where r\ represents the average of rj + 5r] under the assumption (Srj) = 0. The density operator and 
its entropy must depend on the size and type of the noise. If {5r] 2 ) <C {if) 2 the essential effects are 
mostly due to mutual mixing between a few closest eigenstates whose strength for a given eigenstate 
z strongly depends on rj. The places of abrupt structural changes (quantum phase transitions) can be 
easily identified as pronounced maxima of von Neumann entropy [53] . 
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The noise-induced density operator Pi{rf) can also be associated with an equivalent thermally pop- 
ulated system (canonical density operator), which makes it possible to introduce an analog of specific 



heat Ci = Tr[pjln 2 pj] — Sf defined for each state % at given r\ [53]. For a linear hamiltonian (23) and 
small-amplitude noise one obtains an approximation 



^)-«lnW)El™f. <69) 



which, as seen from Eq. (28), represents a direct measure of the parameter- induced variation of wave 
functions. Indeed, for the IBM ground state the "specific heat" introduced in this way yields very 
similar values as the "specific heat" defined through the first derivative of the wave-function entropy 
Wo {5 \v) in the U(5) basis [62]. 

As follows e.g. from Eqs. (28) and (69), the effects of mixing for a given state i are strongly 
enhanced in a vicinity of avoided crossings of level i with neighboring levels. In these cases, the 
energy denominator that appears in the perturbative expansion of the hamiltonian eigenstates makes 



the perturbation V extremely efficient. Recall that hamiltonian (23) can be written in the obvious 
self-similar form H(rj + 57]) = H(r]) + SrjV for an arbitrary point r\. The crossings are avoided for 
levels with the same symmetry quantum numbers because such states generally yield nonzero matrix 
elements (ip^V^j). This prevents a statistically meaningful occurrence of exact degeneracies in spectra 
(corresponding to fixed symmetry quantum numbers) driven by a single (real) parameter rj, although 
actual crossings may occur as exceptions with a vanishing statistical weight (or more frequently in 
integrable systems with hidden quantum numbers). In any case, avoided or unavoided crossings of 
low-energy levels including the ground state seem to hold the key to understanding of both first-order 
and continuous QPT phenomena [52,75]. 

To describe the influence of level crossings near quantum critical points in systems with increasing 



size, it turned out useful to extend the domain of the control parameter rj in Eq. (23 ) from real to complex 
values [82, 108, 134, 137]. The reason for this is the fact that — unlike the real-?] case — the hamiltonian 
eigenvalues do cross in the complex-?] plane. If the complex crossing occurs close to the real axis, one 
observes a sharp avoided crossing in the real-?/ level dynamics. Convergence of the complex crossing 
to real r\ with an increasing size of the system implies that at the corresponding place the evolution of 
energies and wave functions becomes nonanalytic. 

Since the complex-extended hamiltonian H(rj) is no more hermitian (its eigenvalues are complex), 
the points rjk satisfying Ejjrjk) = Ej{r]k) may be called nonhermitian degeneracies. An n-dimensional 



hamiltonian of the type (23) has in total \n{n — 1) complex conjugate pairs of such points. They are 
either diabolic points (a nongeneric type of crossing in complex ??) [176] or branch, so called exceptional 
points (the generic type of degeneracy) [134,177,178]. In a diabolic point, the two eigenvalue sheets 
just touch each other with a conical topology. The branch point, on the other hand, connects the 
two Riemann sheets in a more tricky way, forming a singularity which can be locally described as the 
complex square root. The system of n Riemann sheets associated with solutions of the eigenvalue 
equation becomes entangled. Only the branch points located on the Riemann sheet that for real r\ 
corresponds to the lowest energy eigenvalue can influence the ground-state phase transitional behavior. 
This, however, represents a severe problem because to numerically assign individual branch points to 
their respective Riemann sheets is a very difficult task even for moderate dimensions. 

A method to bypass this problem has been developed and tested in the IBM phase transitions by 
Cejnar, Heinze, and Dobes [82,108]. Its basic idea is that branch points located on the selected Riemann 
sheet near the real parameter axis can be detected indirectly, using only the real energy eigenvalues 
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Figure 5: The evolution of the peak of the ground-state "specific heat" from Eq. (70) at the IBM second-order 



QPT with increasing boson number N. The log-log plots show the peak height (left), the peak inverse width at 
half maximum (middle), and the product of the peak height and width (right). The vertical scale in all panels 
is irrelevant. The linear dependence in the rightmost panel represents a power- law decrease of the peak area, 
indicating a vanishing value of the "latent heat" Q. Adapted from Ref. [108]. 



Ei{rj) for Inxq = 0. To this end, we introduce a quantity 



CM 



n— 1 d-q 2 



Y,^\EM ~ Ei{v)\ 



(70) 



!(.&) 



where n is the dimension of the relevant subspace, counting only the states with the same symmetry 
quantum numbers (for the ground-state considerations, the relevant subspace typically consists of states 
with zero angular momentum, but there may be also other conserved quantum numbers to be taken 
into account, e.g. seniority). As shown in Ref. [108], the quantity Ciirf) measures the proximity of 
branch points located on the Riemann sheet of the ith state to the real parameter axis at i]. In fact, 
it can be derived from a two-dimensional electrostatic analogy, associating the ith sheet branch points 
with charges in the complex plane. Ci(rj) then represents the first derivative of the repulsive Coulomb 
force acting at the place r\ of the real axis: C oc -z^F. If the zth sheet branch points get close to the 
real axis, the change of force at the corresponding place is large and Ci(rf) will show a peak. This peak 
becomes singular in the QPT critical points, where the distance of branch points from the real axis 
asymptotically vanishes. 



At the same time, the quantity (70) can be shown to play a similar role in quantum phase transitions 

In the first-order transition, the QPT peak of 
For a continuous transition, the limiting 



as the specific heat in thermal phase transitions [82, 108] 
Ci(rj) converges to a 5- function located at the critical point r\ c 
behavior of Ci(rj), although still singular, yields a vanishing integral over an infinitesimal vicinity of 
the critical point: Q = lim^o j nc _ e C dr] = 0. This is illustrated in Fig. 5 which shows the boson 
number dependence of the peak height (panel a) and the inverse width at half maximum (panel b) of 
the "specific heat" Co (77) for the IBM ground state at the critical point of the second-order QPT between 
0(6) and U(5) limits. One observes that the peak width decreases faster than the height increases [cf. 
panel (c)], which leads to the estimate of Q = in the N —> 00 limit. This is in agreement with the 
continuous character of the transition studied. Note that for thermal phase transitions, Q would be 
nothing but the latent heat of the transition. 

In thermodynamics, the role of branch points is played by zeros of the partition function in com- 
plex extended temperature. As discussed by Yang and Lee in their 1952 seminal papers [179,180], 
complex zeros of the partition function are crucial for fundamental understanding of the thermal phase 
transitions. The criteria derived for thermal phase transitions of various orders in terms of complex 
zeros [181] are essentially the same as those for QPTs and branch points in the complex parameter 
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plane. Therefore, the above outlined considerations represent a powerful analogy between thermal and 
quantum phase transitions. 

4.3 Signatures involving excited states 

Although the definition of quantum phase transitions is related to the ground state, it is clear that 
there must also exist some QPT induced effects in the behavior of energies and transition amplitudes 
for low-lying excited states. In IBM, such spectroscopic observables can be approached via so called 
quadrupole shape invariants [57,166], which represent expressions based on nth-order moments (n > 2) 
of the quadrupole operator Q x = Q in the ground state: 

92 = (ot\(Q • QM) , qs = (ol\[QQQ] {0) \oi), g 3 = (ot\(Q ■ Q)(Q ■ QM) , ■ ■ ■ (71) 

The mean field expectation values for these invariants can be written in terms of quadrupole shape 
variables (3 and 7, but at the same time, equivalent expressions can be derived with the aid of sums 
over products of reduced matrix elements (0^||Q||2^~) for E2 transitions between individual + and 2 + 
states [166]. Since only the states with small phonon numbers contribute significantly to these sums, 
the shape invariants represent an important tool for identifying suitable QPT signatures in the spectra 
of low lying levels. 

Among the other spectroscopic signatures of the IBM shape transition across the extended Casten 
triangle we name the following: (a) The ratio R4/2 = E(4*)/ ' E{2f) between excitation energies of the 
first 4 + and 2 + states. It takes the harmonic vibrator value R4/2 ~ 2 near the U(5) limit, the rotor 
value ~ 3.33 near the SU(3) and SU(3) limits, and the value ~ 2.5 near the 0(6) limit, (b) The 
quadrupole moment Q(2f) oc (2^||Q||2^) of the first 2 + state, which in the mean field approximation 
is proportional to the ground-state equilibrium value of j3. This measure is related to the expectation 
value (O^ln^lO^) = (rid)o oc f3 2 of the d boson number operator in the ground state [proportional to 
B(E2, 2f — > 0^)], which is however insensitive to the sign of /3 and therefore not suited for the description 
of prolate-oblate transitions, (c) The strength .B(E2,2 2 f — > 2]*~) for the transition from the second to 
the first 2 + state, which is close to zero near the SU(3) and SU(3) limits and nonzero in the spherical 
and 7-unstable phases, (d) The strength p 2 (E0, 0^ — ► Of) of the EO transition between the first excited 
+ and the ground state, which is also related to the <i-boson content of the ground state [75]. (e) The 
ratio £4/2 = B(E2,4f -»• 2^)/ J B(E2,2 ] ^ -> Of) of transition strengths between the lowest L = 0,2,4 
states, which exhibits a peak close to the critical point [78]. (f) Two-particle transfer intensities for 
the ground-state to ground-state transition and transitions from the ground state to excited + states, 
which show characteristic patterns in the critical region [109]. 

An effort has been spent to identify observables suitable for a direct distinction between first- and 
second-order QPTs in finite- N numerical data. Such observables would turn out important in numerical 
studies of systems with a phase transition of unknown nature. In the IBM framework, the quantity 

^2 = M^MO^ -(0+M0+)] (72) 

was proposed as such a measure by Iachello and Zamfir [69] . Indeed, while the above discussed ground 
state average ^(n^) = v\ rises from v\ ~ (spherical phase) to v\ > (deformed phase) with no 
particular difference between the first- and second-order transition, 1/2 is much more sensitive to the 
transition type. For the first-order transition it exhibits an abrupt, sign changing crossover from 1/2 > 
(spherical phase) to z/ 2 ~ (deformed phase), whereas for the second-order transition the change is 
smoother. The key for understanding these observations is the fact that the first-order quantum phase 
transition can be locally interpreted as a sharp anticrossing of the two lowest energy eigenstates with 
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the same symmetry quantum numbers [58]. Indeed, this is how the wiggling behavior of v 2 in the 
first-order transition comes about. The quantities v\ and v 2 can be related to the expectation values 
(r 2 )^, of the squared radius (real space) in the states involved [33], particularly u 2 oc <V 2 ) + — (^ 2 ) + 
measures the isomer shift between the first two + states. In some cases, v 2 has to be replaced by 
v 2 = jj [(^i\ n d\^f) ~ (Qf\ n d\0t)] ; which is proportional to the experimentally available isomer shift 
S(r 2 ) = (r 2 ) 2 + — (r 2 ) +. Note that the isomer shifts are closely related to the expectation values of 
the effective EO transition operator in the states involved, so that there is a link to the EO strengths 
discussed above. 

In a very recent work by Bonatsos, McCutchan, Casten, and Casperson, a new quantity was proposed 
to distinguish the first- and second-order QPT in the IBM framework. It is the energy ratio E(6i)/(0 2 ) 
between the first 6 + and second + states, which takes value 1.5 in the U(5) limit, ~ in the SU(3) 
limit, and 1 in the 0(6) limit. A comparison of some other spectroscopic observables for first- and 
second-order transitions has been presented by Pan, Draayer, Luo, and Zhang [64, 87] and by Zhang, 
Hou, and Liu [112]. 

4.4 Symmetry related approaches 

Essential contributions to the interpretation of low energy spectroscopic data at and around the quantum 
critical points in nuclei have been brought by novel generalizations of the dynamical symmetry concept. 
One of these approaches attempts to ascribe specific dynamical symmetries to the first- and second-order 
critical points. This might seem surprising as the critical point itself is just in between two competing 
symmetries, so one would expect rather disorder than a symmetry at this place. Nevertheless, it turned 
out that these so called critical point symmetries are very successful in the description of spectroscopic 
data. The first specimen of this type was discovered within the fermion dynamical symmetry model 
[141]. The S0(7) dynamical symmetry of that model is located at the phase transition between different 
ground state shapes, similarly as in the IBM the 0(6) dynamical symmetry constitutes the critical point 
between prolate and oblate shapes when proceeding along the SU(3)-SU(3) transition [54]. 

A systematic, though approximate approach to the construction of the critical point symmetries was 
presented by Iachello in the framework of the geometric model [146-150]. Hereby use was made of a 
strongly simplified Bohr hamiltonian which for the spherical-prolate transition led to an approximate 
critical point symmetry denoted as X(5) [147], and for the spherical to 7-unstable transition to the 
symmetry called E(5) [146]. Both these symmetries focused on the /3 degree of freedom, for which an 
infinite square well potential was imposed, and were based on an approximate separation of (3 and 7 
dynamical modes. As a result of these simplifications, virtually parameter free solutions were obtained, 
which prompted an extensive search for experimental examples. 

Since the critical point descriptions are embedded in the framework of the geometric model, they do 
not reflect the finite number of nucleons in real nuclei. Finite size effects have been treated in the IBM 
framework by Leviatan and Ginocchio for the E(5) case (second-order transition) [61] and by Leviatan 
for the X(5) case (first-order transition) [84]. In both cases, the method employed was based on the 
IBM intrinsic states projected onto the appropriate symmetry quantum numbers and on the variation 
of the corresponding energy functional after the projection. In the X(5) case, an additional two level 
mixing calculation had to be carried out as a consequence of the phase coexistence in the first-order 
transition. An effective deformation that results from these calculations constitutes the structure of 
the excited bands at the critical points. However, as explicitly demonstrated by Arias and Garcia- 
Ramos et al. [63,85] for the second-order transition, the critical point symmetries based on the original 
approximation of the square well potential in j3 have only a limited validity in the IBM framework even 
for an infinite number of bosons. This is so because the second-order critical point potential in IBM is 
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essentially a pure quartic oscillator, V oc [3 A . We will therefore not review the X(5) and E(5) symmetries 
in detail here. The interested reader can consult recent overviews by Caprio and Iachello [150], and by 
Casten and McCutchan [23]. 

Another symmetry related approach to quantum phase transition is based on so called partial dy- 
namical symmetries. These symmetries express situations, when only a certain subset of properties 
corresponding to the full dynamical symmetry is retained by the system, whereas the remaining proper- 
ties are not. In particular, all or a part of states may be characterized by a subset of symmetry quantum 
numbers, or a part of states by all quantum numbers [182-186]. Situations like this are particularly 
relevant in systems with coexisting regular and chaotic features, but they may also capture essential 
features of the phase coexistence in first-order QPTs. 

The application of this concept to the first-order critical point in the IBM relies on the possibility to 
increase the size of potential barrier between the coexisting minima within the phase transitional region. 
The two minima, whose existence is inherent to the first-order transition, had been pointed out as a 
potential source of nuclear isomerism [44] and later discussed in a similar context of coexisting spectral 
structures [48,51]. However, for a typical IBM hamiltonian these minima are too shallow [51,78]. A 
method for constructing an IBM hamiltonian with a decent barrier separating the coexisting phases has 
been proposed by Leviatan [94] and recently used by the same author to demonstrate the relevance of 
partial dynamical symmetries at the critical point [106] . The method makes use of the decomposition 
of the IBM hamiltonian into the intrinsic and collective parts [187]. While the intrinsic part generates 
an energy surface with two minima and can be tuned to achieve an arbitrary barrier, the collective 
part is composed of kinetic terms which do not affect the shape of the energy surface. The dynamics 
at the critical point can be again described by an effective deformation determined by variation after 
projection and a two-level mixing calculation. 

The application of this method to the IBM first-order transition led to the identification of a specific 
partial dynamical symmetry at the critical point [106]. Distinct subsets of solvable states, character- 
ized by different types of good dynamical symmetries, were shown to coexist in the system governed 
by the critical hamiltonian. The two subsets of states can be constructed from the U(5) and SU(3) 
dynamical symmetry chains, respectively, in agreement with the limiting dynamical symmetries of the 
transition and with the intuition that the first-order spherical-deformed transition generates coexisting 
vibrational and rotational spectral structures. In fact, also the second-order phase transition in the 
IBM is characterized by a partial dynamical symmetry, although of a different type than the first-order 
transition [all states retain a part of quantum numbers due to the underlying 0(5) symmetry, as dis- 
cussed above] . While the latter feature is most probably specific to the model, the applicability of the 
partial dynamical symmetry with two coexisting subsets of solvable states to generic first-order QPTs 
constitutes an interesting subject for further research. 

The last symmetry-related concept we want to discuss in this subsection is based on so called 
quasi dynamical symmetries. This term should not be confused with partial dynamical symmetries 
discussed above. In a partial dynamical symmetry, a fixed realization of a given algebra is partly 
imprinted in the spectrum (in one of the above senses). In a quasi dynamical symmetry, in turn, the 
symmetry is completely broken for all states, but for a part of the spectrum this happens in a highly 
organized manner which creates an illusion that the symmetry is preserved. Rowe, who coined this 
term, describes it as follows [71]: "Quasidynamical symmetry is an expression of the possibility that a 
subset of physical data may exhibit all the properties that would result if the system had a symmetry 
which, in fact, it does not have." The mathematical formulation of this concept relies on so called 
embedded representations [188], which can be associated with particular (coherent) linear combinations 
of equivalent irreducible representations of a certain Lie algebra that give rise to a new irreducible 
representation. For systems with quantum phase transitions it often appears that two incompatible 
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Figure 6: Amplitudes for an expansion of the lowest eigenstates with L = 0, 2, 4 (iV = 100) in the SU(3) basis 
just before the phase transition to the spherical phase (only a part of amplitudes is shown). A highly coherent 
pattern of mixing illustrates the SU(3) quasi dynamical symmetry. Adapted from Rosensteel and Rowe [83] . 



quasi dynamical symmetries can be associated with the two competing phases. The quasi dynamical 
symmetry of type I starts breaking only if the system enters the phase transitional region, from where 
it may emerge with a quasi dynamical symmetry of type II. 

A good example is the SU(3) quasi dynamical symmetry of the shell and related models [189,190]. 
The exact Elliot SU(3) dynamical symmetry is badly broken by spin-orbit and major shell mixing in- 
teractions, which leads to a complete fragmentation of actual hamiltonian eigenstates among individual 
SU(3) irreps. Nevertheless, the mixing amplitudes preserve a very high degree of coherence, so that 
observable quantities (transition strengths) retain the basic signatures of the unbroken SU(3). 

A similar effect was shown by Rowe, Turner and Rosensteel [78,83] also for the first-order transitional 
path between the U(5) and SU(3) limits of the IBM. An example of the decomposition of eigenstates 
in the SU(3) basis is shown in Fig. [6] The three lowest yrast states shown in the figure at a value of 
the control parameter just before the phase transition to the spherical phase have essentially the same 
amplitudes in the SU(3) basis. In the SU(3) limit, the yrast states are all classified by the same SU(3) 
quantum numbers (A,/i) = (2N, 0). Although the SU(3) dynamical symmetry is apparently broken 
(the states are widely spread in the symmetry basis), the mixing has a highly coherent form. In fact, 
the yrast states on the deformed side of the phase diagram can be derived from an intrinsic state of 
the same form as in the SU(3) limit. The amplitudes evaluated at only a slightly different value of 
the control parameter — just after the phase transition — exhibit a dramatically different pattern, with 
the SU(3) coherence lost. The quasi dynamical symmetry on the U(5) side of the transition can be 
described in terms of the random phase approximation. 

Even a more comprehensive analysis guided by the idea of quasi dynamical symmetry was presented 
for the IBM second-order transition between U(5) and 0(6) limits [70,71], which has a direct counterpart 
in the geometric model [81]. We will return to these results briefly in Subsec. 7.4 
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5 Comparison with experimental data 

If the theoretical signatures of quantum phase transitions are to be compared with experimental observ- 
ables in real nuclei, there appears an essential difficulty: The parameters that induce realistic nuclear 
shape transitions, i.e. proton and neutron numbers Z and A/", do not vary in a continuous way. One has 
to reconcile with the fact that phase transitional regions cannot be studied in arbitrary detail by fine 
tuning of the control parameters. Whether one observes an ideal transitional specimen nucleus between 
two distinct shapes depends on the particular discretization of the relevant section of the nuclear chart. 
Essentially, shapes of atomic nuclei are determined by the changing shell structures of protons and 
neutrons. The presence of a shell closure drives the nucleus towards a dominance of the spherical 
shape. Although the detailed shell configuration represents the microscopic origin of the shape, often a 
simplification can be done when one deals with nuclei far away from (sub)shell closures. The quantities 
that then mostly determine the equilibrium shape of a given nucleus are the numbers N p and N n of 
valence protons and neutrons, or the corresponding holes if the nucleus is located above the midshell. 
A good empirical criterion for the onset of deformation was proposed by Casten et al. [191,192]. It 
reads as 

p -^fk^' (73) 

where P (derived from "promiscuity") expresses an average number of interactions of each valence 
nucleon (or hole) with nucleons (holes) of the other type. For P > 5 the nucleus can be expected to be 



deformed in its ground state, while for P < 5 it is more likely spherical. The approximate condition (73) 
demarcates expected regions of shape-phase transitions across the nuclear chart, as shown in the left 
panel of Fig. [7] [192]. One notices that the regions with P « 5 close to stability are rather scarce. 



Notably, the most important transitional region is centered around M = 90 and at the Os and W 



nuclei. It turns out that the condition (73) is valid close to the U(5)-SU(3) transition, but it is not 
sufficient when a nucleus is close to the U(5)-0(6) transition: For example, 184 W has P = 5.1 but the 
fit of spectroscopic observables places it near the 0(6)-SU(3) leg of the symmetry triangle. 

Since the focus of this review is centered more on the theory side, we present below only an outline 
of relevant experimental issues. The reader interested in more detailed information can consult recent 
reviews [23,24]. 

5.1 Nuclear masses 

Already in their pioneering work on shape-phase transitions in atomic nuclei, Dieperink, Scholten, and 
Iachello pointed out the importance of nuclear mass measurements for the study of shape transitions [9] . 
As the masses depend on binding energies, changes in the absolute energy of the ground state lead to 
kinks in an otherwise smooth behavior. For even-even nuclei these effects together with abrupt changes 
in nuclear radii are the only direct observables since the ground state alone has no electromagnetic 
moments or transitions. 

A good quantity to track changes in nuclear masses is provided by the extraction of two-nucleon 
separation energies, which are extracted as mass differences 

S 2n {Z,M) = M(Z,Af-2) + 2m n c 2 -M(Z,Af), S 2p {Z,N) = M(Z-2,Af) + 2m p c 2 -M(Z,Af) , (74) 

with M(Z,]\T) the mass of a nucleus with Z protons and M neutrons, and m p , m n the masses of 
nucleons. Theoretically, the separation energies can be predicted from the IBM ground-state energies 
Eq for boson numbers N + 1 and N at values rj + Srj and r\ of the control parameter associated with the 
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Figure 7: Left (adopted from McCutchan et al. [192]): a part of the nuclear chart with values of -R4/2 = 
E(4f)/E(2f) and the P = 5 contour, see Eq. (73). The stable nuclei are shown in gray. Right (adapted 



from Garcfa-Ramos et al. [65]): the two-neutron separation energies of Nd, Sm, Gd and Dy isotopic chains. 
Experimental values (full symbols) are shown together with the IBM fit (open symbols). 



given pair of nuclei. Specifically, 

S 2m = E (N+l,ri+6Ti)-E (N,Ti) 



dN 



E (N,r ] ) + 6r ] -§-E (N,r ] ) 



dif' 



(75) 



where • stands for n or p. Since E depends on N quadratically (due to one- and two-body interactions), 
the first term on the right-hand side yields a linear dependence on N and hence a smooth contribution 
to the separation energy. An irregular contribution is created by the second term: The first-order 
transition, with a discontinuous -§-Eo, leads to a discontinuity of the separation energy. The second- 

order transition, with a discontinuous fr^E , induces only a change of slope of the separation energy 
curve at the critical point. 

As an example, the right-hand panel of Fig. [7] compares the experimental values of the two-neutron 
separation energies in the Nd, Sm, Gd and Dy isotopes to the results of IBM calculations performed in 
Ref. [65]. One clearly observes a "discontinuity" (smoothened by finite-size effects) indicating that at 
M = 90 a first-order spherical-to-deformed phase transition is occurring in each of the isotopic chains. 

A systematic study of phase transitions and their effect on separation energies and masses was 
performed by Garcia- Ramos, De Coster, Fossion, and Heyde in Ref. [55]. These authors find a very 
high sensitivity to the IBM parameters around the second-order phase transition, which may indicate 
that the description of energies and electromagnetic transition rates is not a sufficient criterion in this 
region. Compared to the first-order transition, the effect of a second-order phase transition on separation 
energies is by about an order of magnitude smaller (a change in the slope instead of a discontinuity). 
It is therefore very important that nowadays very accurate mass measurements with AM/M ~ 1CT 7 
become possible on stable as well as exotic nuclei [193] . The access to exotic nuclei is essential to verify 
phase transitional predictions in the P m 5 regions far from the valley of stability. 

5.2 Low-energy collective states 

The shape of the ground state is strongly present in the properties of the low-lying excited states. In 
order to unambiguously determine the shape types from the data, Casten et al. [22,50,191,194] used 
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Figure 8: The excitation energy of the 4f state against one of the 2f state for the nuclei with E(2f) < 600 keV. 
The spherical-deformed transition is observed as the change of the slope between the indicated vibrational and 
rotational values at E(2f) ss 150 keV. 

some universal dependences obtained from spectroscopic observables for a large number of isotopes. 
The best known example is shown in Fig. [8j where the excitation energy of the first 4 + state is plotted 
against the one of the first 2 + state. All data together clearly yield a universal curve which shows two 
different slopes. At low excitation energies the slope is approximately 3.33, which indicates rotational 
motion, while at high energies the slope becomes approximately 2, indicating a spherical vibrator. 
The first-order shape-phase transition occurs where both slopes cross. One can thus deduce a typical 
excitation energy for the first excited state as a signature for the spherical-deformed phase transition. 
Among the other indicators of shape-phase transitions are the electromagnetic properties, especially 
the quadrupole moment of the 2f state and strengths of E2 transitions between excited states (cf. 
Subsec. 4.3). For the prolate-oblate transition, the change of the sign of (3 must induce the change 



of the sign of the 2 { quadrupole moment. A more difficult task is to observe and distinguish in data 
the second-order phase transition, as this relies on the determination of the critical exponent from 
rather scarce data points. In the following we focus on systematic studies of chains of nuclei along the 
U(5)-SU(3), U(5)-0(6), and SU(3)-SU(3) transitional paths. 

5.2.1 Spherical to prolate transition 

Early work on a systematic comparison with experimental data for the transitional region between U(5) 
and SU(3) was done by Scholten, Iachello and Arima [195]. They studied in particular the Sm isotopes 
where the spherical-deformed transition is very well pronounced. These isotopes and in particular 
152 Sm, which is an almost ideal first-order transitional nucleus, formed the basis of the renewed interest 
in shape-phase transitions. This was triggered again when a new study of 152 Sm showed that the 
transition strength ratio B(E2; 2+ — > 0^)/-B(E2; 2+ — > Of) is extremely small, but could be reproduced 
by a simple IBM calculation [196]. It was then shown that the sharp minimum in -B(E2; It — > Of) is 
located at the first-order phase transition, in the region where spherical and deformed phases coexist [48] . 
This interpretation launched the revival of interest in shape-phase transitions, especially after the 
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Nucleus 


148 Nd 


150 Nd 


150g m 


152 Sm 


154 Gd 


156 Dy 


176 0s 


178 0s 


P 


3.75 


4.4 


4.4 


4.8 


5.1 


5.3 


4.5 


4.6 



Table 1: P factors from Eq. (73) for the isotopes situated close to the first-order transition. 




11 

Figure 9: The number of + states observed below 2.5 MeV as a function of r] (see text). The corresponding 
nuclei are from left to right: 162 Dy, 168 Er, 158 Gd, 176 Hf, 154 Gd and 152 Gd. 



introduction of the X(5) critical point symmetry [147] and the subsequent interpretation of 152 Sm as an 
example of a X(5) nucleus [197]. The interpretation of shape-phase transition caused an intense debate 
which is still ongoing (for a review see Ref. [23]). As an example, the moments of inertia in the bands 
corresponding to both phases are almost the same, contrary to the model prediction [51]. 

Further nuclei located at or very close to the first-order transition were studied experimentally. They 
are the M = 90 isotones 150 Nd [198] and 154 Gd [199]. They also follow the X(5) pattern in energies 
and electromagnetic transition rates for the ground state band. The next candidate, 156 Dy, follows still 
the X(5) energy trend but in its electromagnetic decay it already behaves as a rotor [199]. Besides 
the J\f = 90 nuclei, which are situated directly in the region with P = 5 (see Fig [7]), the neutron 
deficient Os nuclei were recently proposed as other candidates. Lifetime measurements indicate that 
indeed 176 ' l7 " 8 0s exhibit the X(5) symmetry and are well described in the IBM as nuclei at the first-order 
critical point [200,201]. 

A systematic comparison of rare earth nuclei was done by Garcia- Ramos, Arias, Barea, and Frank 
[65]. Chains of Nd, Sm, Gd and Dy isotopes were fitted using the simple consistent-Q hamiltonian as 
well as a complete IBM-1 hamiltonian. Good descriptions for the excitation energies and B(E2) values 
were obtained together with a good description of the two- neutron separation energies (see Fig. [7|). 
Based on the fitted parameters the isotope chains could be placed in the separatrix plane studied before 
within the framework of the catastrophe theory [47] . The results showed that none of the Dy isotopes 
is close to the critical point. The closest ones are 148 Nd, 150 Sm and less clearly 152 Gd, followed by 
150 Nd and 152 Sm. The finding that the latter two isotopes are not so close as 148 Nd and 150 Sm is a bit 
surprising in view of the low P-factor for 148 Nd, see Table [lj From the values shown in the table one 
may deduce that the first-order phase transition occurs when P is about 4.5 rather than 5. 



Using the simple IBM hamiltonian of the type (56), McCutchan et al. [192] performed systematic 
fits of Gd, Dy, Er, Yb and Hf isotopes, which allowed the location of these isotopes in the triangle. It 
was found that the phase transition is crossed inside the triangle for 154 Dy and 156 Er, while 152 Gd lies 
on the U(5)-SU(3) leg before the phase transition and 154 Gd lies close to it, but after the transition. 
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The Yb and Hf isotopes do not cross the phase transition on the U(5)-SU(3) side. 

Another experimental observable was proposed by von Brentano et al. [75] who observed that at 
the first-order transition from spherical to deformed shapes a rapid increase of the E0 strength for the 
®t ~~ y Oi" transition takes places. A recent measurement in 154 Sm confirmed this behavior nicely [202]: 
p 2 (E0; 0+ -> Of) = 0.020(2) for 150 Sm, 0.051(5) for 152 Sm, and 0.096(42) 154 Sm. 

Finally, as shown e.g. in Ref. [52], at the first-order phase transition the spectrum of excited states is 
maximally compressed. Using the (p, t) reaction at the highly sensitive Q3D spectrometer in Garching, 
an extended search for + states in several even-even nuclei in the rare earth region was performed [102]. 
It was possible to identify the + states up to very high energies. Fig. 6 shows the numbers of observed 
+ states below 2.5 MeV for the nuclei studied in Ref. [102] against the fitted value [203] of r\ in the 
hamiltonian parametrization (56). One clearly observes an increase of the + state density at the critical 
value 7] = 0.8 (corresponding in the fit to 154 Gd). Note that here no corrections were applied for the 
changing boson number and the fitted \ parameter. 

5.2.2 Spherical to 7-unstable and prolate to oblate transitions 

While there is ample evidence for atomic nuclei situated at the first-order phase transition between 
spherical and prolate deformed shapes, the situation is less clear for the second-order phase transition. 
Despite the fact that the U(5)-0(6) transition in the IBM has been throughoutly studied from the 
theoretical side, only few examples of nuclei located along this leg of the symmetry triangle are known. 
Early work [206] indicated that the neutron-rich Ru and Pd isotopes are well described by this transi- 
tional class. More recently, the Ru isotopes were reanalysed by Frank, Alonso, and Arias [56] to locate 
the second-order critical point. Using two neutron separation energies, level energy systematics, and 
B(E2) values these authors concluded that the chain of Ru isotopes in between 98 Ru and 110 Ru can be 
described as a U(5)-0(6) transition with the critical point situated at 104 Ru. 

After the introduction of the E(5) critical point symmetry [146], the nucleus 134 Ba was proposed 
as a candidate for the second-order critical point nucleus by Casten and Zamfir [204]. In this work, a 
comparison of B(E2) and energy ratios with an IBM calculation at the second-order phase transition 
was made, yielding a very good agreement. A subsequent systematic search in the nuclear database 
indicated that 128 Xe might be another good candidate for the second-order critical nucleus [205]. Very 
recently, the Zn isotopes were studied and 64 Zn was proposed to be located close to the second-order 
transition [207]. 

For the SU(3)-SU(3) first-order transition the critical point takes place at the 0(6) dynamical 



40 



symmetry. Although oblate nuclei, or transitions from 7-soft to oblate shapes, are relatively rare in the 
known part of the nuclear chart, they do occur in the Pt and Hg isotopes close to the 0(6) nucleus 
196 Pt. In Ref. [60], three different signatures were investigated: the 4+ to 2f energy ratio R4/2, the 2f 
quadrupole moment Q(2f), and the transition strength B(E2; 2\ — > 2f). In an attempt to span a larger 



part of the extended Casten triangle, these quantities are plotted in Fig. 10 for nuclei ranging from the 
well deformed prolate rotor nucleus 180 Hf up to 200 Hg. For each nucleus the \ value from the quadrupole 



operator Q x was fitted using the hamiltonian (56) with r\ = and the respective value of N. Clearly 
all experimental observables indicating the prolate-oblate phase transition are observed [68]. The very 
rapid change of the quadrupole moment might explain why 196 Pt has a nonvanishing quadrupole moment 
while still being a good 0(6) nucleus. Surprisingly, the Hg isotopes do not resemble vibrational or shell 
model like structures whom have the R4/2 ratio around or below 2. 

6 Quantum phase transitions in IBM extensions 

As mentioned in Sec. [3j the interacting boson model with just two types of bosons represents the 
simplest member of the IBM family. Extensions are possible in several directions: The distinction 
of proton-proton and neutron-neutron pairs leads to the two-component IBM-2, while mixed proton- 
neutron pairs are included in more sophisticated versions of the model (IBM-3 and -4). The description 
of odd nuclei can be achieved when valence states of a single fermion are added into the Hilbert 
space within the IBFM. To describe other types of excitations, one may consider non-scalar bosons of 
different angular momenta (p, f,g,...) together with (or instead of) the d boson; this is also relevant if 
the IBM approach is applied in molecular or other types of physics. The inclusion of higher order boson 
interactions, external rotation, or different particle-hole configurations represent some other possibilities 
to refine the model. Each of these possibilities implies particular modifications of the model phase 
structure. Some of these modifications will be outlined in this section, although the research in this 
field is still in progress. The study of phase transitions in the IBM extensions is interesting not only 
from the viewpoint of concrete applications, but also because these models provide case examples of 
the influence of additional degrees of freedom on the QPT phenomena. 

6.1 Proton-neutron degrees of freedom (IBM-2) 

In the version of the interacting boson model discussed up to now (IBM-1), no distinction was made 
between proton-proton and neutron-neutron pairs. In order to make the model more realistic, it is 
essential to introduce bosons of the proton and neutron type. It is done in the two-component version 
of the model, the IBM-2 [33]. Hereby it is implicitly assumed that the shell model orbits which the 
protons and neutrons occupy are different, as is the case for medium mass and heavy nuclei. The IBM-2 
gives a successful phenomenological description of low-energy collective properties of virtually all such 
(even-even) nuclei. The IBM-3 and 4, which are applicable in lighter nuclei, are not discussed here. 

In the IBM-2, the total number of bosons N is the sum of the neutron and proton boson numbers, 
N v and N n , which are both conserved separately. The spectrum generating algebra of the IBM-2 is a 
product Uy(6)®U 7r (6), consisting of neutron and proton generators b vln} b v ii m > and frj^fe^m/, respectively. 
The model space is the product [N v ] x [N n ] of symmetric representations of U^(6) and U ff (6). The most 
general (N u , N n ) conserving rotationally invariant IBM-2 hamiltonian will not be presented here, as it 
has 21 parameters. Instead, a strongly simplified hamiltonian is used. A comprehensive review of the 
model and its implications for nuclear structure can be found e.g. in Ref. [208]. Here we highlight three 
main features. 
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The first is that the existence of two kinds of bosons offers the possibility to introduce a so called 
F-spin quantum number [209], which is rather similar to isospin. Since F — |, the bosons can be in 
two possible charge states: M F — — | for the neutron type boson and M F = +| for the proton type. 
The F spin is formally defined by the algebraic reduction 

U(12) D U w (6) <g> U(2) 

I I I (76) 

[iV] [JV - /, /] [N - f, f] 

with 2F = N — 2f the difference between the Young tableau labels that characterize both \J un (6) and 
U(2). The algebra U(12) consists of the generators b'i m b p 'i' m i, with p, p' = v or n, which also includes 
operators that change a neutron boson into a proton boson or vice versa. Under the U(12) algebra, all 
bosons behave symmetrically, hence the representations of U V7r (6) and U(2) are identical. 



The SU(2) algebra associated with the F-spin, which is a subalgebra of U(2) in (76), consists of the 
diagonal operator F z = \{N^ — N v ) and the raising and lowering operators F± that connect neutron and 
proton bosons. It is clear that all IBM-2 hamiltonians must commute with F z . The hamiltonians which 
also commute with F 2 have the F-spin as a good quantum number. The classification and analysis of 
the F-spin conserving dynamical symmetries of the IBM-2 show that the three symmetry cases of the 
IBM-1 are recovered. These are the U^(5), CU(6) and SIU(3) limits of the IBM-2 [210]. 

The second important aspect of the model is that it predicts states which are additional to those 



found in IBM-1. Their structure can be understood in terms of the F-spin classification (76) for F- 
symmetric hamiltonians. The states with the maximal F-spin value, F = |iV, are symmetric in U^(6) 
and represent exact analogues of the IBM-1 states. The next class of states has F = |iV — 1, and these 
are no longer symmetric in U I/7r (6) but belong to its representation [N — 1,1]. Such mixed-symmetry 
states were studied theoretically in 1984 by Iachello [211] and since then they have been observed in 
many nuclei (see e.g. Ref. [212]). Of particular relevance are the 1 + states, since these are allowed 
in IBM-2 but not in IBM-1. From the geometrical analysis performed in the limit of large boson 
numbers [213] it emerges that the mixed-symmetry states correspond to linear or angular displacement 
oscillations, in which the neutrons and protons are out of phase, in contrast to the symmetric IBM-2 
states for which such oscillations are in phase. The occurrence of mixed-symmetry states was first 
predicted in the context of geometric two-fluid models in vibrational [214] and deformed [215] nuclei, 
in which they appear as neutron-proton counter oscillations, pictorially referred to as scissors modes. 

The third important feature offered by the IBM-2 is the possibility to get triaxial shapes besides 
the axially symmetric ones. This happens when proton and neutron fluids exhibit different types of 
deformation — prolate and oblate. Hereby a new dynamical symmetry can be constructed, denoted by 
SU* 7r (3), for which the respective x parameters are in absolute value equal to ^- but have opposite 
signs [216,217]. This possibility gives rise to the Dieperink tetrahedron, which has an extra dimension 
compared to the Casten triangle, and to a new, triaxial shape phase of the model. Note that if all 
neutron-proton deformation combinations are taken into account, there are in fact four and not just two 
SU(3) limits: prolate-prolate, SU(3), oblate-oblate, SU(3), prolate-oblate, SU(3)*, and oblate-prolate, 
SU(3) *. Among these combinations, however, only two are inequivalent. 

The analysis of shape-phase transitions in the IBM-2 is much more complex than the one in IBM-1. 
It was done independently and at the same time by Arias, Garcia- Ramos, and Dukelsky [72], and by 
Caprio and Iachello [73,92]. The classical limit of the IBM-2 can be obtained using a two-component 
coherent state, 
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Figure 11: The phase structure of the IBM-2 hamiltonian (78) in the space of parameters (, xs = \{Xtt + Xv), 
and xv = \{Xtt ~ Xv)- The surfaces represent first- and second-order phase transitions between spherical and 
deformed, and axially symmetric and triaxial shapes. Adapted from Caprio and Iachello [92]. 



which is used to construct the energy surface for a given hamiltonian as a function of parameters 
.>,,,! ., Those can then be interpreted as 4 deformation parameters jj n , r y in &v l lv l and 2x3 



anm and a vm 



Euler angles. In contrast to the IBM-1, the Euler angles do partly affect the energy surface through 
the relative orientation {01,62,63) of the proton and neutron ellipsoids. Therefore, the IBM-2 energy 
functional needs to be minimized in 7 parameters. 

As mentioned before, a tractable study of phase transitions using the complete IBM-2 hamiltonian 
(with 21 parameters) is hardly feasible. Hence only some simplified hamiltonians, resembling those used 
in the IBM-1, were studied. The following hamiltonian was used in both Refs. [72,73] 



F*-*"(C) = (1 - ()Un d * + n dv ) + (^ [-{Q*- + Q*>) ■ (Q*« + Q*>)] 



(7S 



with N the total boson number. The form (78) is a direct generalization of the IBM-1 hamiltonian (56). 



It contains (apart from an overall scaling factor that we do not account for) 3 control parameters (, Xm 
and Xv- The hamiltonian belongs to the type for which it was shown [213] that the global minimum 
only occurs for aligned proton and neutron deformations, reducing the number of order parameters from 
7 to 4, i.e. the /3 and 7 shape variables of both fluids. 

The phase structure of hamiltonian (78) can be studied analytically for £ = 1, where the O w (6), 
SUj,,r(3), SU* 7r (3), and the remaining two equivalent SU(3) dynamical symmetries are located. The 
ground state at O y7r (6) represents a 7-unstable deformed configuration, at SU V7r (3) a 7-rigid deformed 
configuration, prolate or oblate for both fluids, whereas at 811*^(3) the ground state is triaxial since 
both fluids have opposite types of deformation (with the symmetry axes orthogonal to each other) [213]. 
A second-order phase transition between axially symmetric and triaxial deformed shapes occurs at the 
boundary determined by the following equations [73] , 



9f/M-l)+A,(2/^-l)(^ + l) 
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N w 
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with f3 p expressing the deformation parameter of the neutron (p = v) and proton (p = n) ellipsoids. The 
location is strongly dependent on the ratio N^/Ny except at the 0^(6) limit which is always located 
at the second-order transition between 7-soft and triaxial deformation. 

A numerical analysis of the full parameter space showed that the above phase separatrix extends from 
( = 1 to ( < 1, where it forms a surface of second-order phase transitions between axially symmetric and 
triaxial deformed phases. This plane is terminated at the spherical-deformed phase transition, which 
also forms a surface close to the 11^(5) dynamical symmetry. The spherical-deformed phase transition 
is of the first order except at the intersection with the axisymmetric-triaxial phase transition, which 
happens at \ v = Xn — 0. The whole phase diagram in the symmetry tetrahedron is depicted in Fig. [TTJ 
A very detailed account, including an analysis of more complex hamiltonians can be found in Ref. [92]. 

A closely related problem was recently analyzed by Iachello and Perez-Bernal [116] in the context of 
molecular physics, within the so called vibron model for coupled benders [218]. The model represents 
a generalization of the 2D vibron model, which is a lower-dimensional analog of the sd-IBM with the d 
boson replaced by a two-component boson r [219]. The spectrum generating algebra of the 2D vibron 
model coincides with U(3), while the extension investigated in Ref. [116] concerns coupled systems (such 
as molecules of the ABBA type) with the spectrum generating algebra U(3)®U(3). A simple ground- 
state phase diagram of three basic molecular configurations was drawn in the plane of the hamiltonian 
control parameters. 

6.2 Systems with bosons and fermions (IBFM) 

One particularly important extension of the simple interacting boson model concerns odd-mass nuclei. 
It is achieved by considering, in addition to the s and d bosons, a fermion coupled to the core with an 
appropriate boson-fermion interaction. The resulting interacting boson-fermion model (IBFM) [35] is 
thus a specific version of the particle-core coupling model, which has been widely used in nuclear physics 
to describe odd-mass nuclei [165]. The characteristic feature of the IBFM is that it lends itself very well 
to a study based on symmetry considerations whereby certain classes of boson-fermion hamiltonians 
can be solved analytically. The basic building blocks of the IBFM are N bosons, in terms of which the 
even-even core states are modeled, and fermions occupying a set of single-particle orbits with various 
angular momenta j. Low- lying collective states of an odd-mass nucleus with 2N + 1 valence nucleons 
are approximated as a single fermion coupled to iV-boson states. A particularly attractive feature is the 
similarity in the description of even-even and odd-mass nuclei. This has given rise to the development 
of a supersymmetric model [220] . 

The IBFM hamiltonian contains bosonic and fermionic parts and a boson-fermion interaction: 
-£/bf = H-q + Hp + Vbf- The bosonic hamiltonian can be taken directly from the IBM. Because only one 
fermion is coupled to the boson core, the fermion hamiltonian can be written as Hp = ^ . ejUj, where 
rij is the number of fermions in orbit j and 6j its single-particle energy. The boson-fermion interaction 
is assumed to have a two-body character: 

Vbf = J2 v Wf [( & ? x 4) (J) x & x a /) (J) ] f ■ ( 80 ) 

iji'fj 

This form, however, is too general for phenomenological applications. A simplification is given by: 
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Here : • : indicates normal ordering and Qb denotes the boson quadrupole operator, Eq. (44). The 



exchange interaction Vg F takes account of the fact that the IBM bosons have an internal structure, and 
this leads to the exchange effects. The coefficients k can be related to the occupation probabilities and 
quasi particle energies derived from a BCS calculation (see [35]). 

The spectrum generating Lie algebra associated with the IBFM is [Ub(6) ®Uf(^)]. The usual boson 
algebra Ub(6) describes the collective core excitations, while Uf(^), with Q = J2j(^J + 1) denoting the 
size of the fermionic single-particle space, corresponds to the fermion degrees of freedom. The algebra 
contains two sets of generators: in terms of bosons &} m &z' m ', and in terms of fermions aj m a,,-/ m / (both 
types written most favorably in the angular momentum coupled form). The existence of analytically 
solvable IBFM hamiltonians relies on isomorphisms between some boson and fermion algebras. A simple 
example is provided by the angular momentum algebras defined separately for bosons and for fermions. 
The former consists of the angular momentum operators Lb which generate the boson algebra Ob (3) and 



occurs in the lattice (46). The fermion angular momentum operators Jp are generators of an equivalent 
Lie algebra denoted as Slip (2). An essential difference between these algebras is that while the bosons 
couple only to integer angular momenta, the fermions can form also states with half-integer momenta. 
The summed operators Jbf — Lb + Jf generate the boson-fermion algebra SUbf(2) of the total angular 
momentum, which can be integer or half-integer depending on whether the number of fermions is even 
or odd. Since Ob (3) has to be a subalgebra of Ub(6) and SUp(2) one of Up(^), it follows that SUbf(2) 
is a subalgebra of the IBFM spectrum generating algebra. The required form of all IBFM dynamical 
symmetries is the following: 

[U B (6) ® U P (fi)] D • • • D SU B f(2) . (82) 

The subscript BF is usually omitted and the total angular momentum algebra is denoted as Spin(3). A 
comprehensive review of various dynamical symmetries of the IBFM can be found in Ref. [35]. Numerous 
experimental examples are given in Ref. [222]. 

It is clear that in absence of fermions the IBFM reduces to the standard IBM of even-even nuclei. 
This naturally leads to a super ■symmetric extension of the IBFM that allows a simultaneous description 
of even-even and odd-mass nuclei. Let us stress that a necessary condition for such an approach to 
be phenomenologically successful is that the energy scales for bosonic and fermionic excitations are 
comparable. Nuclear supersymmetry treats spectral properties of different nuclei as arising from a 
single boson-fermion hamiltonian and a single set of transition operators. This unification is achieved 
by embedding the IBFM spectrum generating algebra into a superalgebra U(6/f2): 

U(6/fi) D [U B (6) <8> Up(fi)] . 

I ' i i (83) 

[AO [N] [1 M ] 

Here, the supersymmetric representation [J\f} of U(6/0) imposes symmetry in the bosons and anti- 
symmetry in the fermions, and contains the [U B (6)®U F (f2)] representations [N] x [1 M ] with M = N+M. 
Thus, a single supersymmetric representation contains states in even-even (M = 0) as well as odd-mass 
(M = 1) nuclei. Nuclear supersymmetry was originally postulated by Iachello [220] as a symmetry 
among pairs of nuclei. Subsequently, it was extended to quartets of nuclei, where the neighboring 
odd-odd nuclei are also incorporated [223]. 

To probe the ground-state phase transitions in the IBFM, we consider below a fermion in orbits 
with angular momenta j = |, |, and |. The relevant dynamical algebra is [Ub(6) <8> Uf(12)] and an 
isomorphism between the boson and fermion algebras is established by introducing pseudo-spin and the 
pseudo-orbital angular momenta for the fermion: Sp = \ and L F = 0, 2, respectively. This leads to 
a reduction of the spectrum generating algebra to [Ub(6) <8> Uf(6) <g> Uf(2)], which allows to combine 
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Figure 12: Absolute energies of the + states in even-even nuclei (a) and of the corresponding ^ states in 
odd- .A nuclei (b) along the prolate-oblate shape transition. The calculation was done using the hamiltonian 
(84) with N = 10, r\ = 0, and \ variable. Adapted from [77]. 



both U(6) algebras and to construct the three dynamical symmetries of the standard IBM [224]. The 



embedding (83) results in the U(6/12) supersymmetric scheme. While this is sufficient for a theoretical 



study, the practical application is limited to those nuclei where the above fermionic shells are important, 
e.g. the A =70-80, the Rh-Ag and the Pt-Os-Hg mass regions. 

The first study of quantum phase transitions in odd-,4 nuclei was due to Jolie, Heinze, Van Isacker, 
and Casten [77]. It relied on a supersymmetric extension of the results obtained for even-even nuclei. 



The simplified hamiltonian (56) was generalized to 



tf*^) oc (1 - t?U (-Q* F 



Q£p)+»A[Ubp(5)] 



34) 



with Ci[Ubf(5)] the first-order Casimir operator of Ubf(5) and Q^ F the quadrupole operator of the 



form (44) but using generators of Ubf(6). The use of U(6/12) supersymmetric scheme allowed to 



extend the Casten triangle to odd-mass nuclei. In Ref. [77], the application of Eq. (84) was discussed 
to the prolate-oblate phase transition on the SUbf(3)-Obf(6)-SUbf(3) line (77 = 0). 



Figure 12 compares the absolute energies of the + states of the N = 10, M = system to the ones 
of Jbf = \ states in the N = 10,M = 1 system. In the even-even system (M = 0), the prolate-oblate 
phase transition is associated with the kink of the lowest energy at 0(6) and the sequence of energies 
is symmetric around 0(6), although the wave functions are different [54]. A similar behavior is found 
for the odd-A. case (M = 1). Nevertheless, a closer inspection reveals some interesting differences. 
In the even-even case, except at the dynamical symmetries with x = or ±^, all states undergo 
avoided level crossings indicating the absence of conserved quantum numbers (besides the total angular 
momentum). This is not so in the odd-A. calculation where one observes real crossings of several states. 
The origin of these crossings follows from the conservation of the angular momentum Lbf — Lb + Lp by 



hamiltonian (84) which results from the coupling of the total boson angular momentum to the pseudo- 
orbital angular momentum of the fermion. Lbf is conserved in classifications contained in the lattices 
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Figure 13: Comparison between experimental and theoretical negative-parity states of odd-neutron isotopes 
between 191 Os and 199 Hg. Adopted from [77]. 
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Whatever embeddings are taken in the dotted part of (85), Lbf remains a good quantum number and 
the same holds also for the pseudo-spin Sp. A detailed analysis shows that the essential ingredient is 
the proper choice of the single-particle energies of the pseudo-spin doublet j 



I and | [77]. 



In Ref. [77], the hamiltonian (84) was applied to odd-^4 nuclei between 191 0s and 199 Hg as a logical 



extension of the fit done for the even-even cores (see Fig. 10 ) [68] . It was shown that the hamiltonian 
needs to be extended with C 2 [Ubf(6)] and C 2 [Spin(3)], which do not affect the phase transitional 



behavior but allow a fine tuning of energies. The result is shown in Fig. 13 Although the hamiltonian 
is very simple, and although the parameters were kept constant across the region (after having adopted 
the values for r\ and x from the fit to the even-even partner), the overall structural changes, in particular 
the decreasing trend and compression of the low-lying levels with decreasing mass, were reproduced 
rather well. 

Critical behaviors of odd-^4. nuclei have been recently studied for the Ubf(5)-Obf(6) second-order 
transition by Alonso, Arias, Vitturi, and Fortunato [86,101,107,110]. These analyses were motivated 
by the fermion critical point symmetry E(5/4) introduced by Iachello [149] as an extension to odd 
masses of the E(5) critical point symmetry [146]. The E(5/4) symmetry is obtained in the geometric 
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model when coupling a j — | fermion to the E(5) core, but a similar system can be obtained also in 
the IBFM [86] . While the applicability of critical point symmetries is restricted only to critical points 
of phase transitions, the IBFM makes it possible to describe the whole transition from vibrational to 



7- unstable nuclei. In Ref. [86], a specific form of the interaction (80) was used with only the quadrupole 



term in the boson-fermion interaction (81). Imposing \ — for the boson quadrupole operator in both 
the boson hamiltonian and the boson-fermion interaction, the total hamiltonian conserves in addition to 
the 0(5) quantum number also the Spin(5) quantum number along the whole transitional region. From 
the comparison between the E(5/4) and IBFM approaches at the critical point a moderate agreement 
was obtained. 

In a subsequent work [101], the boson-fermion interaction was extended with more terms conserv- 
ing the Spin(5) symmetry. In addition to the electromagnetic transitions also single-particle transfer 
intensities were calculated. However, even the extended hamiltonian did not remove the discrepancies 
between both approaches. The differences can be traced back to the square well form of the potential 
of the geometric model assumed in critical point symmetries, in contrast to the /3 4 dependence of the 
critical IBM potential [85]. 

Because the fermion space with only a j — | orbit is rather restricted, Alonso, Arias, and Vitturi 
have extended recently both models to allow the fermion to occupy also orbits with j = |, |, and 
I [107,110]. For the critical point symmetry this leads to the E(5/12) structure and for the IBFM 



it allows to use the U(6/12) supersymmetric approach as discussed before. The hamiltonian (84) was 
used to study the second-order phase transition in odd- A nuclei, posing now as a constraint x — 0. 
The variable parameter 7] controls the transition from Ubf(5), at r\ = 1, to Obf(6), at rj = 0. The 
energy spectrum and electromagnetic transition rates were compared with the critical point symmetry 
E(5/12) [110]. In this case, a better agreement between both approaches is obtained and signatures of 
the critical behavior are found in energies and B(E2) ratios of low- lying states. However, no comparison 
with experimental data is made in these works. 

An extension of the coherent state formalism to odd-mass systems with a single fermion orbit was 
recently presented by Liu [111]. Under the neglect of the 7 degree of freedom, the energy functional was 
derived and studied for a hamiltonian with only a quadrupole term in the boson-fermion interaction. 
In this case, a somewhat surprising dependence of the positions of critical points rj c and Xc on the 



interaction strengths k from Eq. (81) was obtained. Further studies including a proper treatment of 
the 7 degree of freedom are needed to verify these results. 

The study of shape-phase transitions in odd-^4 nuclei is more complicated than the one in even- 
even nuclei, also because the level densities are much higher in the former case. To investigate mixed 
systems of bosons and fermions is nevertheless of great interest as an increasing number of Bose-Fermi 
mesoscopic systems becomes available for experimental studies. 

6.3 Other extensions 

6.3.1 Configuration mixing (IBM-CM) 

Besides the shape-phase transitions atomic nuclei may also exhibit shape coexistence. This phenomenon, 
invoking multiple shapes present in the low-energy spectrum, was predicted and observed in many 
spherical nuclei near magic shells. Particular configurations observed in these nuclei can be linked to 
the occupation of specific up- and/or down-sloping orbitals, coined "intruder orbitals" , which allowed 
for a simple understanding of the shape coexistence in the works of Heyde, Wood et al. [225,226]. The 
method put forward by these authors has been used to predict shape coexistence e.g. in the Sn nuclei 
around mass number A = 116, and in the light Pb nuclei below A = 196. 
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The difference between shape coexistence and shape-phase transitions is not always very clear, 
especially from an experimental point of view [74]. In a shape-phase transition, the number of basis 
states is preserved through the transition as progressing from one limiting case into the other. The 
eigenstates of one limit are spread out among the actual, intermediate eigenstates, which finally become 
the eigenstates of the other limit. At any instant, the eigenstates form a complete basis of the same 
Hilbert space. On the other hand, typical situations implying shape coexistence appear when a set of 
extra states (outside of the model space that is regularly considered for low-lying configurations) drops 
in excitation energy and produces new (intruder) low-lying states. The Hilbert space is extended as 
H = H nor © H int , where the two subspaces correspond to normal and intruder states. Any hamiltonian 
acting in this space decomposes into three parts: the one of normal states, the one of intruder states, 
and the one describing their interaction. In practical IBM applications, this division is implemented 
by coupling the sd boson subspaces with total numbers of bosons N and N + 2. These are used to 
describe shape coexistence between normal (valence) excitations and intruder excitations that include 
2p-2h configurations. Both sets of states can interact with each other via a mixing hamiltonian that 
does not conserve the boson number. 

Both these seemingly different pictures were connected in the approach recently proposed by Frank, 
Van Isacker, and Iachello [96] . These authors have introduced a method for studying phase transitions 
in systems with mixed configurations, which they call type-II phase transitions. As explained above, a 
configuration mixing (CM) hamiltonian consists of parts corresponding to normal and intruder states 
(these form diagonal blocks of the hamiltonian) and their interaction (which appears off the diagonal). 



To cope this structure with the QPT paradigm expressed by the hamiltonian in Eq. (24), a general 
phase transitional CM hamiltonian was considered in the form: 



H{fi,u,A) 



1 - 77)i?i uW 
uW rjH 2 + A 



36) 



Here, Hi and H 2 are operators acting in the space of normal and intruder states, respectively, while W 
represents the interaction hamiltonian scaled by the mixing strength oj. The normal and intruder parts 
of the hamiltonian are driven by an external parameter r\ between the rj = and rj — 1 limits, in which 
only one of the parts is active. An offset A energetically separates intruder and normal configurations. 
To determine the phase structure of the CM system, the concept of an eigenpotential [227] has been 



employed. Hamiltonian (86) is evaluated in the appropriate group coherent states of the respective 
subspaces. These, in the IBM case, may be associated with condensate states of N and N + 2 bosons 
(for Hi and H 2 , respectively). The interaction W couples both types of coherent states, so one obtains a 
2x2 hamiltonian matrix expressed in terms of the coherent state parameters that capture the geometry 
of the problem. The lower eigenvalue of this matrix (the lower eigenpotential) represents the coherent 
state expectation of the ground state energy. The minimization of this eigenvalue in the coherent state 
parameters yields the ground state energy as a function of the hamiltonian control parameters, which 
can be subject to the phase transitional analysis. 

In Ref . [96] , an example has been analyzed when the normal and intruder configurations are described 
by the U(5) and 0(6) dynamical symmetries [228]. The coherent state averages of the hamiltonian 



blocks in Eq. (g took forms (H x ) = p 2 /(l + p 2 ), (H 2 ) = J[(l - P 2 )/(l + P 2 )] 2 , and (W) = to (since 
W = s* s' + S •aH-h.c). If uj — A = 0, the hamiltonian (86) describes two noninteracting configurations 
with a pair of degenerate (3 = and (5 > minima, and the entire line r\ G [0, 1] is critical. For uj > 0, 
the phase structure becomes more complicated and strongly dependent on the offset A. If A is small, 
the first-order critical curve in the plane rj x uj is surrounded by spinodal and antispinodal lines that 
converge to a single curve with increasing uj. For uj — > oo, the transition becomes of the second order; 
this limit is equivalent to the normal U(5)-0(6) case. For larger values of A, however, a region with 
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a single spherical minimum develops near 77 = 1 which keeps growing until it eventually extincts the 
phase coexistence regions. For A > | the ground state eigenpotential has only one minimum, either 
spherical or deformed. 

It should be stressed that these structures are essentially different from those obtained within the 
Landau approach to phase transitions. This is because the above procedure, involving the diagonal- 
ization of the potential matrix, does not lead to the ground state energy in the form of a polynomial 
in the order parameter. General consequences and other examples of this generalization are still to be 
investigated. A characteristic signature of the type-II first-order phase transition seems to be a large 
phase coexistence region. In the original paper [96], the Sr and Zr nuclei at M = 58 were proposed 
as candidates for this type of transition. The above approach was recently applied by Hellemans, Van 
Isacker, De Baerdemacker, and Heyde [113] to the case of so called U(5)— Q x ■ Q x mixing in nuclei, 
which includes both U(5)-0(6) and U(5)-SU(3) types of phase coexistence. Let us note that similar 
models can be applied also in molecular physics. 

6.3.2 External rotation (IBM cranked) 

Interesting results have been obtained by analyzing the phase transitional properties of the system of 
IBM bosons in a rotating frame. The so called cranking approach, which imposes an external rotation 
onto an ensemble of interacting particles, is a standard tool to describe moments of inertia and other 
collective aspects of rotational motions. The cranking analysis of nuclear quadrupole shapes performed 
on the fermionic level by Alhassid et al. [21, 229, 230] in the framework of the Landau theory led to 
the determination of a shape-phase diagram for hot rotating nuclei. The cranking approach within the 
IBM was pioneered by Schaaser and Brink [231,232] and by Dukelsky et. al. [233,234]. The mean field 
phase transitional analysis of the cranked IBM was performed by Cejnar [59,62]. 



The procedure is similar to that discussed in Subsecs. 2.1 and |3.2[ It consist of evaluating the 



expectation value of the cranked hamiltonian H' = H — lu-L (where a; is a vector of the cranking angular 



frequency and L the angular momentum operator) in a condensate state (51 ). An important difference, 
however, concerns the form of the condensate: its coefficients cannot satisfy the usual condition a +m = 
(— ) m a*_ ml resulting from the hermicity of the Cartesian matrix atij, since in such a case the contribution 
of the cranking term uj ■ L vanishes. Instead, one has to assume some modified conditions [232] that for 
a concrete choice of rotation along the x axis can be represented by the following parametrization [62] : 

«o = /9cos7, a±i = 4=/? sin 7 sin <5 , a± 2 = 4s/3 sin 7 cos 5 . (87) 

The generalized condensate states for uj 7^ then carry angular momentum while the standard results 
for uj = are recovered by setting 6 = 0. One has to be careful since (3 and 7 do not any more carry the 
standard shape interpretation. The geometry involved in the generalized coherent state must be read 
out from the expectation value of the quadrupole operator in these states. Minimization of the cranked 



energy functional with respect to /3, 7, and 5 for the simple IBM hamiltonian (56) leads to the phase 



diagram depicted in Fig. 14 , Note that the finiteness of angular momentum within the IBM (L < 2N) 
puts an upper bound on the physical domain of the cranking frequency uj. In the U(5) limit, e.g., the 
limiting frequency coincides with the boundary of the spherical phase, see Ref. [62]. 

The effect of external rotation may be related to the effect of an external magnetic field. Indeed, 
Fig. [14] resembles a magnetic phase diagram of a superconductor. The same type of phase diagram is 
expected to hold for superfluid and normal phases in nuclei [18, 19], but in the context of quadrupole 
shape phases it represents a modification of the earlier shell model based result of Ref. [21]. In particular, 
the present approach predicts a certain critical rotational frequency that induces a transition from the 
spherical to deformed phase (similarly as there exists a critical magnetic field for a transition from the 
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Figure 14: The phase diagram of the cranked hamiltonian (56) in the plane rj x uj. Cranking frequency u is 



expressed relative to O from the energy scaling factor a = Ml. The upper frequency limit is indicated by the 
dashed line. Adapted from [59]. 



superconducting to normal phase), which is in contrast to the previous result. This becomes interesting 
in connection with the recent proposal of Regan et al. [235] to interpret experimental data on the yrast 
backbending effect in spherical nuclei as signatures of a rotation induced spherical-to-deformed shape 
phase transition. A discussion of experimental yrast band data in 108_114 Cd, 104 ~ 108 pd, 100 ' 102 Ru, and 
100 Mo within the cranked IBM framework can be found in Ref. [76]. 



6.3.3 Other types of bosons and higher-order interactions 

Other IBM extensions, whose phase structure has so far been explicitly analyzed only to a limited 
extent, include in particular (a) the incorporation of bosons of other types (these can either replace, or 
supplement the existing d boson) and/or (b) the addition of higher-order interactions (to supplement 
the present one- and two-body interactions). 

The case (a) with the d boson replaced by another type of boson, let us denote it as b, is well studied 
for lower dimensions. In fact, the Lipkin model corresponds to b = t (a scalar or pseudoscalar boson). 
Such modifications of the standard IBM-1 are experimentally relevant in molecular physics, where b = p 
(a boson with spin / = 1, hence with 3 angular momentum projections) corresponds to the 3D vibron 
model [37,38], while b = r (a two-component boson with no angular momentum content) to the 2D 
vibron model [218,219]. Related approaches have been also applied in condensate state physics [42,236]. 
The phase structure of such models in general has been considered by Cejnar and Iachello [103]. It 
turns out that all two-level boson models of this type, with U(n) as a spectrum generating algebra, 
exhibit a second-order QPT between "spherical" and "deformed" phases, characterized by dynamical 
symmetries U(n — 1) and O(n) (the spherical phase can be unambiguously characterized by the ground 
state having the form of s-boson condensate). The first-order transitions, however, are a much scarcer 
spice. If the angular momentum conservation is required [for models with b corresponding to a boson 
with an even number of components we consider rotations generated by the group 0(2) instead of 
0(3)], the occurrence of a first-order transition of the above type is limited to spectrum generating 
group dimensions n = 2, 6, 10, . . . (corresponding to cases when boson b has even angular momenta). 

Phase transitions in the 2D vibron model (n = 3) were discussed by Pan et al. [88]. A detailed study 



was recently presented by Perez-Bernal and Iachello [115]. As already mentioned in Subsec. 6.1, the 
latter calculations were even extended to a two-component version of the model (analogous to the IBM- 
2) [116]. Phase transitions in the 3D vibron model (n = 4) were investigated by Van Roosmalen [45] 
and by Arias, Dukelsky, Dusuel, Garcia- Ramos, and Vidal [89,104], using techniques beyond the mean 
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field (cf. Subsec. 4.1). The case n = 5 was briefly discussed in Ref. [103], but so far this model has 
no application. From the QPT viewpoint, all these results are essentially equivalent to those of the 
standard IBM (n = 6) between U(5) and 0(6) dynamical symmetries (except that the IBM yields a 
richer system of excited states). 

The problem becomes increasingly difficult starting from the cases b = / and especially b = g, i.e. 
bosons with angular momenta I = 3 and 4. For I = 4, both first- and second-order transitions between 
spherical and deformed phases are supposed to take place. To determine the phase structure seems 
to be complicated already on the mean field level since the energy functional depends in general on 
2/ + 1 variables (radius (3 plus 2/ hyperspherical angles) while the elimination of Euler angles (which 
is essential in the I = 2 case) is not feasible for I > 2. Moreover, in nuclear physics the bosons 
with higher angular momenta are useful mostly as supplements of the dominant quadrupole degrees of 
freedom, which are carried by the system of s and d boson. One should therefore deal with situations 
in which p, f,g, . . . bosons are considered together with s and d. The phase structure then depends on 
a very large number of parameters and severe simplifications are necessary. Although analytical and 
phenomenological studies of additional non-quadrupole degrees of freedom within the IBM exist (see 
e.g. the coherent state analysis in Ref. [237] and the spdf calculations in Ref. [238]), no comprehensive 
QPT analysis is available for such systems so far. 

The situation is similar also for the case (b) of the above classification: the higher-order interactions. 
Inclusion of such interactions within otherwise standard sd-lBM-l may lead to non-trivial effects (cf. 
Ref. [185]). The most general energy functional can be fabricated just from various powers of expressions 
j3 2 and j3 3 cos 37 [since for I = 2 there exist only two independent scalar couplings of shape coordinates, 



see Eq. (54)]. For instance, the energy functional up to three-body terms looks as follows 

A0 2 + B(3 3 cos 3 7 + C{3 4 + D{3 5 cos 3 7 + E(3 6 cos 2 3 7 + F(3 6 
(A7) ~ (1 + /? 2 ) 3 ' 

where the coefficients A, . . . , F depend on the concrete hamiltonian. 



The form (88) generates a new type of phase transition, namely the one connected with the angle 
variable 7, as recently discussed by Iachello [148]. To make it clear, let us rewrite the energy functional 
in the form V(@, 7) = f((3) — g((3) cos 37 + h(/3) cos 2 37, where the functions / and g can be easily read 



off from Eq. (88), and denote the point of global minimum as (j3 m , j m ). For h((3 m ) = the minimum 



in angle variable is apparently at 7 m = [for g((3 m ) > 0] or at 7 m = |7T [for g((3 m ) < 0], both cases 
corresponding to axially symmetric (deformed) shapes. An increase of h(/3 m ) above the critical value 
— \g{@-m) eventually turns 7 = (or |tt) into a saddle point, creating a new minimum < 7 m < |7T. 
This corresponds to a second-order phase transition from axially-symmetric to triaxial shapes [148]. 
Such a situation was discussed for a concrete choice of three-body IBM-1 hamiltonian by Jolos [79,80]. 
A special type of a three-body hamiltonian was considered by Thiamova and Cejnar [97]. It was 
chosen to preserve the 0(6) dynamical symmetry in the whole parameter space, although the underlying 
0(5) symmetry could be broken [185]. The three-body interaction generated by the term [QoQoQo}^°\ 
where Q represents the 0(6) quadrupole operator, induces the occurrence of the /3 3 cos 37 term in 
the energy functional and hence a crossover from 7- soft to 7-rigid shapes. There is no shape-phase 
transition in angle variable in this model (since the term oc cos 2 37 is absent), but the phase structure 



is very much reminiscent of the extended Casten triangle of hamiltonian (56) with its three phases and 
a triple point in between. 

Let us conclude this section by noting that the general phase structure of k-body IBM hamiltonians 
was not yet determined, and this not even for k = 3. 
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7 Excited-state quantum phase transitions 

So far we have been probing the system of interacting bosons at about zero temperature. However, 
phase transitions driven by interaction parameters can also appear at finite temperatures. In the last 
part of this review, we turn our attention to possible extensions of the QPT description to states with 
an arbitrarily high degree of excitation. In contrast to the traditional thermodynamical treatment, 
we will mostly deal with individual states instead of ensembles. We will show that in the infinite 
size limit excited states throughout the whole spectrum may exhibit nonanalytic evolutions with the 
control parameter at a certain sequence of "critical" points. For this phenomenon we use the name 
"excited state quantum phase transition" (ESQPT) [100, 117]. Concrete investigations in this part 
are mostly restricted to the IBM between U(5) and 0(6) dynamical symmetry limits, but we believe 
that continuation of work in the directions outlined here may bring new results also in more general 
situations. 

7.1 Finite temperatures 

We start with some remarks on the thermodynamical extension of the ground-state QPT. A system in- 
teracting with a heat bath at fixed temperature T is naturally described [1] by the Helmholtz free energy 
F = (E) — TS, where (E) = Tr(Hg) is the energy average and S = Tr(— glng) the entropy associated 
with the quantum density operator g (we set the Boltzmann constant to unity, so the temperature is 
measured in units of energy and the entropy is dimensionless). Maximizing the entropy while keeping the 
energy average fixed leads to the familiar thermal population of states p%(T) = ^4^ exp(— T~ l E i ) 1 hence 
q = Qo(T) = Y(f) ex P(~T~ 1 H), with Z(T) = ^2 i exp(—T~ 1 E i ) = Trexp(— T~ X H) being the canonical 
partition function. The equilibrium value of the free energy is obtained by substituting g = go(T) in 
the above formula, which yields 

F (T) = (E) T - TS (T) = -T In Z(T) . (89) 

The following well-known relations can be derived for the derivatives of the equilibrium free energy with 
respect to temperature, 



A. 

c)T 



F (T) = -S (T) , ^F (T) = -£ [(E 2 } T - {E%] , (90) 



T 2 C(T) 

where C(T) is the specific heat. Since the energy dispersion in the second equation is never negative, 
the equilibrium entropy Sq(T) is a nondecreasing function of temperature. The thermodynamical phase 
transition of nth order (n = 1, 2, . . . ) corresponds to a situations when the value Sq(T) itself (for n — 1) 
or its (n— l)th derivative (for n > 1) jumps discontinuously (in the thermodynamical limit) at a certain 
transitional temperature T c . If the (n — l)th derivative (n > 1) gets infinite at T c , the phase transition 
is called just continuous (with no order in the Ehrenfest sense). 

At zero temperature, the equilibrium free energy coincides with the ground-state energy: F (0) = 
(E) = E . The ground-state QPT, expressed through the dependence of E on control parameters, may 
therefore be a special case of a more general phenomenon that resides in the space of the model control 
parameters and temperature. In infinite lattice models with finite range of interactions, the extension 
of QPTs observed at T = to T > is of general interest and numerous results can be found in the 
literature [11,12]. In finite many-body models, results on T > phase transitions are not abundant. 
This is in spite of the fact that an analytic approach to describe thermodynamical phase transitions in 
pseudospin models as well as in some non-spin algebraic models was formulated by Gilmore [7, 8] . 
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Basic steps of Gilmore procedure can be briefly outlined using the hamiltonian based on the U(2) 



spectrum generating algebra, i.e. the Lipkin model of Subsec. 2A We may use the representation in 
terms of a spin-| array (counting N constituents) with infinite-range interactions. First, one makes 
clear that in the mean field approximation (which for the present class of models becomes exact in the 
N — > oo limit) the thermal iV-body density operator can be expressed as a product of 1-body density 
operators, i.e. in the form g( N > = Qi®Q2<S>. . ■ Qn = (f?) 0JV 5 where Qh — Q are identical operators acting in 
the 2-dimensional Hilbert space of fcth particle. Each g^ can be written in the usual form g — |l + a ■ a, 
with 1 standing for the unity matrix and (a x , a y ,a z ) for Pauli matrices. Polar coordinates in the 
decomposition (a x ,a y ,a z ) = (r sin 6 sin (ft, r sin 9 cos 0, r cos 6) represent a convenient set of parameters 
to determine the state g( N \ The minimization of the free energy F^ N ' = Tr[Hg( N '] — TTr[— g^ In g^] 
with respect to these parameters yields the equilibrium value F for each set of external parameters 
and temperature. The free energy for large values of iV can be written as F^ N ' = Nf(r,9,(f)). For 
T = the minimization gives r = |; then g( N > represents a pure state and the method reduces to the 
condensate state procedure used for the ground state. For T > 0, however, one gets r < |. Phase 

transitions (if any) can be detected as nonanalytic evolutions of _F = Nfo with the model control 
parameters. In this way, a finite-T counterpart of the ground-state phase transition in the Lipkin model 
was detected [7,10]. 



With a hamiltonian satisfying the linear ansatz (23) the equilibrium free energy F depends just 
on two parameters: T and rj. As may be anticipated from the Gilmore method, for finite models the 
description of phase transitions in the rj x T plane will not be essentially different from the Landau 
theory of thermodynamical phase transitions in the p x T plane, where p stands for another intensive 
thermodynamical variable like pressure. One therefore expects that phase transitions will form some 
curves (or parts of curves) in rj x T. Consider a point P on such a curve and at this point mark the 
normal and tangent directions n and t, respectively. For the classification of the phase transition at P, 
discontinuous partial derivatives of Fq along the normal direction n are substantial (partial derivatives 
along t are continuous). However, since partial derivatives along an arbitrary direction can be expressed 
through partial derivatives of the respective order along n and t, the same classification is likely to be 
found if probing the transition along any of the two parameters rj and T. 

Let us explicitly evaluate the first two partial derivatives of F with respect to t] and T. The 
derivatives ^Fq and ^Fo at constant 77 obviously yield Eq. (90). For the derivatives containing t] one 
obtains, 

IF ( V ,T) = (V) T , $F (j,,T) = -HiV~ (V)t]\ + (|U) t , (91) 

where (X)t = 'Yl l iPi{T)X i with X{ = (\l/j|X|\l/j) represents the thermal average of a quantity X. An 
analogous expression can be derived for ^g-F (r],T). Note that the T = limit of these equations 



coincides with Eq. (29), discussed above in connection with the ground state. 



The first equations of both pairs (90) and (91) show that a jump of entropy (in the vertical, i.e. 



temperature direction) is likely to be connected with a jump of the interaction average {V)t (in the 
horizontal, i.e. interaction parameter direction). These are equivalent expressions of a generic first-order 
phase transition at such places in rj x T where the phase separating curve is not parallel with one of the 



axes. Similarly, as seen from the second equations in both pairs (90) and (91), a nonanalyticity of the 



specific heat is expected to have an equivalent counterpart in the expression of the thermal dispersion 
of V and/or the thermal average of the level curvatures jj^iEi- 

In view of the above discussion, one may tend to anticipate that any QPT at 77 = r] c and T = 
smoothly continues to the T > domain, where it can be detected as a thermodynamical phase 
transition at T = T c {rf). Such a statement is the content of a so called crossover theorem, formulated by 
Gilmore for the specific systems he analyzed [7,10]. It should be stressed, however, that this rule cannot 
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be applied in general. The simplest counterexample arises if symmetry constrains disable repulsion of 
some levels and the ground state undergoes an unavoided crossing. Then a T = phase transition 
(typically a first-order one) takes place at the crossing point, with no continuation to T > 0. 

To illustrate this mechanism on a trivial example, let us take a specific IBM hamiltonian with the 
U(5) dynamical symmetry: 

H{rj) = (1 - r])n d + V n s = V N + (1 - 2 V )n d =► V(m /3) = r) + (1 - 2rf) -?—= . (92) 

1 + p z 

Here, n s and n d stand for the numbers of s and d bosons, and V represents the classical potential 



derived from the hamiltonian. The eigenstates of (92) for any value of the control parameter r\ can be 
classified by the U(5) quantum numbers and the spectrum becomes fully degenerate at T] = T] c — |. 
For rj < rj c the energy of the ground state, which is an s-boson condensate, linearly increases with r\. 
At T) — T) c the ground state flips ito the (i-boson condensate and its energy starts to linearly decrease. 
From the viewpoint of the above nomenclature, such a situation would be classified as the first-order 
ground-state quantum phase transition; this being so for an arbitrary (even finite) number of bosons N. 
However, since all levels vary with the control parameter in an analytic (linear) way, the introduction 
of a nonzero temperature makes the free energy an analytic function of both rj and T. The phase- 
transitional behavior detected at T = disappears for an arbitrarily small temperature T > 0. Other 
examples of finite- N ground-state crossings on the 0(6)-U(5) side of the Casten triangle can be found 
inRefs. [66,98]. 

In spite of its success in pseudospin models, the approach based on canonical thermodynamics 
has some drawbacks in general finite systems. In particular, we know that quantal spectra of such 
systems (for finite numbers of particles) have finite numbers of states, so the level density in the upper 
part of the spectrum is a decreasing function of energy. Under such circumstances, the introduction 
of temperature can work as in usual thermodynamical systems only at some low temperatures. For 
higher temperatures, the thermal averaging smoothenes all properties over the whole spectrum and 
disables a direct access to states with higher energies. (The only way to bypass this problem would be 
to consider negative absolute temperatures as in [52].) For this reason, it is usually more convenient 
to use directly the excitation energy instead of temperature. Quantum phase transitions for excited 
states then show up as singularities of the level density and also as nonanalyticities in the evolution of 
individual level energies and wave functions. Since the logarithm of the level density is proportional to 
the microcanonical entropy, the present approach is related to the microcanonical thermodynamics, see 
e.g. Ref. [5]. As canonical and microcanonical results should converge in the proper thermodynamical 
limit, the present approach does not imply any loss of thermodynamical information. 

7.2 Excited level dynamics 



It is clear that evolution of the spectrum of the general linear hamiltonian (23) with parameter r\ is fully 
determined by a "snapshot" containing all energies Ei = Ei(rj) (with i — 0, . . . ,n — 1) and all matrix 
elements Vij = {^(^Vtyjirj)) & t an arbitrary single value of r\. This viewpoint was elaborated by 



Pechukas and Yukawa [129,130], who considered the level dynamics induced by hamiltonian (23) with 



variable r\ as a classical problem of determination one- dimensional motions of n interacting particles. 

A 

drf 



In this approach, parameter rj plays the role of time, while Et are positions and 4-Ei = Ei velocities of 



individual particles. The dynamical equations 

A F - t/ _^L/?-9\^ iJJd Am. - ST v ^Vkj sr^ VikVkj (Qq . 

dr) *i-Vii, dv ^i-*2^ E E > dr, v v- 2^ Ei-Ek ^ Ei-E k [ ' 

j&i) ° k{+i) k(&) 3 
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Figure 15: The grid of IBM-1 states with L = in the plane nj versus energy per boson for N 
(a) and (b) correspond to 0(6) and SU(3) limits, respectively. 



40. Panels 



(which immediately result from ordinary perturbation theory) are analogous to those describing dy- 
namics of a gas of particles interacting through two-dimensional Coulomb force Fij oc qiqj/(xi — Xj). 
The essential difference, however, is the fact that the product of charges g^- is replaced by \Vij\ 2 , which 
cannot be factorized and is also subject to variation due to the third equation. The knowledge of E iy 
Ei and V^- at some 77 = 77^ allows the determination of level dynamics for any 77. 

From the computational viewpoint, there is certainly no advantage in using the system of coupled 



differential equations (93) instead of diagonalizing the hamiltonian (23) for each particular value of 77 



However, Pechukas- Yukawa theory is fruitful in gaining deeper insight into the mechanisms leading to 
QPT behaviors. It shows that whether a given model will exhibit a QPT or not is entirely determined by 
the "initial conditions", i.e. by values of Ei, Ei and V^- at an arbitrary 77 ini . In case when the two phases 
of the model coincide with some dynamical symmetries, the initial point is naturally put to one of these 
symmetries, which is therefore made responsible for the phase-transitional behavior at some i] c 7^ 77^ . 
Moreover, the evolution of quantum energies and interaction matrix elements with 77 is pictured as an 
intuitively imaginable process a la classical mechanics. This makes it easier to guess which types of 
initial conditions are "QPT friendly" and which are not. In particular, a plausible condition for a QPT 
to occur is a fast initial compression of the whole spectrum, which leads to quickly increasing level 
repulsion and may eventually produce a sharp scattering of levels. In the following we show that this 
type of mechanism indeed applies in the IBM case. 



Let us have a closer look on the initial conditions for the IBM simplified hamiltonian (56). A natural 
0, which coincides with the 0(6) dynamical symmetry for \ 



initial point is r\ xx 

or SU(3) dynamical symmetries for x — ± 



2 ' 



0, with the SU(3) 
or is intermediate between these cases for other values of 



X € [— 2 ' + 2 ]• ^ s n °ti ce d in Ref. [98], a large fraction of the spectrum with L = in the 0(6) and 
SU(3), or SU(3), limits exhibits a fast initial compression, having almost a sharp point of focus on the 
U(5) side of the transition. This is true on the whole 77 = side of the extended Casten triangle. 



From the first equation in (93) we can easily derive the condition for a subset S of eigenstates of 
hamiltonian (23) at 77 = r/ in i pointing to a single focus (r}i oc , E{ oc ); it reads as 



(tpi(r] ini )\H(rj ioc )\ipi(r} ini )) = E ioc for i e S 



(94) 
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In our case, with 7/ ini = and the focal point at (rji oc ,Ef oc ) = (1,/), the condition (94) can be written 
as (ipi(0)\rid\ipi(0)) = f N. The value of {nd)i is easy to calculate for the ground state with large N, 
using the coherent state relation (n^)o = /3m/(l + Pm)- I n the 0(6) limit we have (n d ) = |AT, while 

In the 0(6) limit, the condition {na)i ~ ^A" is valid for 



in the SU(3) and SU(3) limits (n d } = \N 

a great majority of excited states with r = and L 

values of seniority and angular momentum, see Fig 



15 a 



0, but it is increasingly violated for growing 



hand, the condition (n^i ~ \N holds for a large portion of the whole L = spectrum, see Fig 



In the SU(3) and SU(3) limits, on the other 

Tsfb). 

This means that the initial compression in the latter case is faster and involves more states, which is 
consistent with the fact that the subsequent QPT is of the first order, i.e. sharper that the one observed 
in the 0(6)-U(5) transition. 



An interesting type of solution of the system in Eq. (93) which may apply under certain initial 



conditions is a formation of a wave propagating through the spectrum. A spectacular example of such 
a phenomenon is shown in panel (a) of Fig. 16 It displays the spectrum of the IBM states with zero 



spin and seniority between the 0(6) and U(5) dynamical symmetries [hamiltonian (56) with x — 0]. 
The "shock wave", which is initiated in the uppermost part of the spectrum at the 0(6) side of the 
transition, keeps moving due to the initial compression of the spectrum and eventually reaches the 
ground state at the QPT critical point r) c = I. This effect was for the first time investigated by Heinze, 
Cejnar, Jolie, and Macek in Ref. [98]. Examples from a broader class of models of the same type have 
been subsequently analyzed in detail by Caprio, Cejnar, and Iachello [117]. 



As discussed in Subsec. |7.4[ the curvatures of individual level trajectories as they pass through the 
wave in Fig. Il67a), i.e. the deceleration and acceleration due to the repulsive force in (93), becomes 



infinite in the asymptotic size regime. Therefore, the wave can be viewed as a manifestation of quantum 
phase transitions for individual excited states (ESQPT) at a running value of the control parameter. 
This interpretation is supported by the analysis of the behavior of branch points in complex extended 



control parameter (cf. Subsec. 4.2). Indeed, the branch points corresponding to individual avoided 
crossings of levels along the wave propagation get closer to the real r\ axis as the number of bosons 



increases. Panel (b) of Fig. 16 shows the "specific heat" from Eq. (70), which measures the proximity of 



branch points on the zth level Riemann sheet to the real axis for a given N. The moving peak coincides 
with the trajectory of the wave in panel (a). The analysis by Cejnar, Heinze, and Macek [108] showed 
that the behavior of the peak for very large A" is consistent with the assumption of a continuous ESQPT 
(the "latent heat" Q is zero). This conclusion is confirmed by semiclassical arguments [100,117] that 



will be outlined in Subsec. 7.4 Thermodynamical consequences may be anticipated from Eq. (91) and 



are probably in agreement with the early analysis by Gilmore [7] (cf. Subsec. 7.1). 



The transition from SU(3) or SU(3) to U(5) does not seem to show any ESQPT effect, at least not 
at such degree of coherence as in the previous case. The exceptional behavior of the [0(6)-U(5)]dO(5) 
transition is probably connected with the fact that the corresponding hamiltonian is integrable but, 
at the same, creates repulsion of levels with the same quantum numbers L and r. Note that the 
level repulsion is not a generic feature of integrable systems since the underlying integrals of motions 
[unlike the 0(5) Casimir invariant in the present case] usually depend on the control parameter [239]. 
The "laminar flow" of levels (due to integrability) and the level repulsion [due to persistent 0(5) 
symmetry] constitute ideal conditions for the formation of a coherent wave. In nonintegrable transitional 
regimes, the chaotic ("turbulent") flow of levels (with numerous collisions of neighboring levels) obscures 
the observation of any collective aspect in level dynamics. So far, the existence ESQPTs was not 
systematically investigated in these regimes. 
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Figure 16: Panel (a) shows a "shock wave" formed in the N = 100 spectrum of IBM states with L = r = 
along the transition from 0(6) to U(5), described in the parametrization (56). The energy is expressed per 



boson. Panel (b) shows the "specific heat" from Eq. (70) evaluated across the spectrum, with x denoting the 
relative excitation ratio (x = for the ground state and 1 for the highest state). 



7.3 Monodromy 

It turns out that the sequence of excited state quantum phase transitions along the 0(6)-U(5) path is 
closely related to so called monodromy. In classical mechanics, monodromy is a topological property 
of some integrable systems that usually indicates the presence of an anomalous (pinched) torus of 
trajectories in the phase space [240]. This anomaly prevents a full analytic description of the system 
(in spite of its integrability) and creates a dichotomic distinction for the types of classical orbits. The 
pinched torus shows up as a point defect in the lattice of quantum states, in which the states are 
organized with respect to their quantum numbers. If a closed loop around the defect is followed, 
the form of an elementary cell in the lattice varies such that after the return to the initial point the 
transformed cell does not fit with the one taken at start. Hence the term monodromy, i.e. "once 
around" . 

The simplest systems that exhibit the monodromy are (a) the spherical pendulum [240,241] and (b) 
particle moving in the sombrero potential [242] . In both these cases, the pinched torus — see Fig. [l7|a)— 
is formed by orbits passing the points of unstable equilibrium (at the north pole of the pendulum sphere 
or at the top of the sombrero potential) with just the critical energy E c at which to reach that point 
takes infinite time. Since this happens at zero angular momentum J z , the corresponding motions can 
be reduced to an effectively two-dimensional phase space, whose intersection with the pinched torus 
defines a classical separatrix. The dynamics on both sides of the separatrix shows qualitatively different 
features: while orbits below the critical energy (inner side of the separatrix) are essentially vibrations 
around the stable equilibrium, those above the critical energy (outer side of the separatrix) extend over 
a larger part of the phase space. In the quantum case, if all angular momenta J z = mh are considered, 
the states below and above the critical energy form different types of patterns in the plane E x J 2 , as 
can be seen — for the sombrero potential — in Fig. [i~7|fo). 

Monodromy in the two-dimensional pendulum and sombrero systems affects bending vibrations of 
some triatomic molecules, such as H2O [242], where the concept was also confirmed on the experimen- 
tal level [243]. The algebraic description of these molecules is possible within the 2D vibron model, 
which is a two-level bosonic model with a scalar boson and a two-component boson [115,218,219], 



cf. Subsec. 6.3.3 A detailed study (including the analysis of excited state quantum phase transitions, 
monodromy, scaling laws and thermodynamical quantities) was recently presented by Perez-Bernal and 
Iachello [115]. The phase transitional properties of the 2D vibron and related models are very close 
to those of the sd-IBM along its 0(6)-U(5) transition. It turns out that in absence of 0(3) rotations 
(L = 0) the seniority quantum number r can be regarded as an analog of the 0(2) angular momentum 
quantum number m considered in the above 2D examples. 



Recall that if the IBM hamiltonian (47) contains no admixture of the SU(3) invariant, the cubic term 



in the classical potential (55) is missing: B = 0. In this case, the critical point at A = separates the 



sombrero and oscillator types of potential (located on the deformed and spherical sides, respectively). 



The full classical hamiltonian in the parametrization (|56|) with x — is given by 

hm = ^ ot + fi + (1 - b) p 2 ] 



:i+W + ( ? 



(5r) - A)(3 2 + r//3 4 

T+J 2 



+ ' n : o 2 , 2 - . (95) 



— v v-~ 

Tvib V 



where T rot and T vib stand for the rotational and vibrational kinetic energies; pp G [0, 00] and p 1 G [0, 1] 
are momenta associated with the shape coordinates (3 and 7, respectively. The radius /3 G [0, 00] can 



be compressed to J3 G [0, \/2] through the relation (3 = (3/y2 — (3 2 [121], while the radial momentum is 
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Figure 17: Panel (a): a schematic view of the pinched torus in the phase space x x y x p r of the sombrero 
potential. Panel (b): the lattice of quantum states in the sombrero potential (adapted from Child [242]). Full 
(broken) lines join points with the same radial (principal) quantum numbers; the circle indicates a crystal 
defect associated with the monodromy. Panels (c)-(e): quantum lattices of IBM states with L = for the 
transition between the 0(6) (77 = 0) and U(5) (77 = 1) limits for N = 40 (adapted from Ref. [99]). 
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transformed accordingly: pp i— > pp G [0, v2]. Eq. (95) then takes the form 

2" 



njjj) = T rot + 2 + ( i_^ 



*K? 



+ (1^-2)^ + (l-r7)/3 4 



(96) 



which is more convenient in the present case since the transformed potential is no longer distorted 
by the /3-dependent denominator. If considering motions with L — 0, the classical phase space is 
generated solely by the vibrational degrees of freedom, described either by polar coordinates (3 and 7, 
or by the corresponding Cartesian coordinates x = /3cos7 and y = /3shi7. The angular momentum 
p 1 = xp y — yp x is then related to the classical limit of the 0(5) Casimir invariant for L = through 



^C 2 [0(5)] d =p 



As shown by Cejnar, Macek, Jolie, Heinze, and Dobes [99,100], the sombrero potential from Eq. (96) 
for r] < r] c = I generates both classical and quantum signatures of monodromy. It was recognized that 
the energy E c = corresponding to the pinched torus separates two different dynamical regimes that 
can be linked, on both quantum and classical levels, to the respective limiting dynamical symmetries 
0(6) and U(5). In particular, the following observations were made: (a) Arrangements of states in 
the Ext quantum lattice below and above E c resemble those associated with dynamical symmetry 
limits 0(6) and U(5), respectively [99]. In the 0(6) type of lattice states form well distinguished chains 
with an approximately quadratic dependence on r, while in the U(5) type of lattice the dependence on 
t is essentially linear, seen Fig. flWc)-(e). (b) The form of wave functions changes from 0(6)-like to 
U(5)-like at E c [117]. The crossover is most clearly seen for r = states. Below E c the average (rid)i 
gradually decreases with energy and the spread of the n^ distribution increases, while above E c these 
trends are reverted, (c) Classical trajectories change their typical forms at E c [99]. This can shown by 
expressing the ratio Ry/p = Ty/Tp of periods associated with 7 and (3 degrees of freedom. Trajectories 
below E c exhibit a broad distribution of Ry/p with the lower bound at ~ 3 and a strong preference for 
large values. This corresponds to flower-like orbits in the sombrero potential. The trajectories above 
E c show a narrow distribution concentrated at values < 3, which with an increasing energy converges 
to the lower limit Ry/p = 2. This corresponds to the bouncing-ball orbits traversing across the central 
maximum. 



7.4 Semiclassical considerations 

The relation of monodromy to nonanalytic evolution of individual r = excited states with the control 
parameter was for the first time noticed [100] in application of the "displaced oscillator" approximation 



of Rowe [70,71] to hamiltonian (56). The approximation is valid on the 0(6) side of the transitional path 
for asymptotic boson numbers, when x = 2^ — 1 can be treated as a continuous variable. Eigenstates 
of the hamiltonian are expressed as conventional wave functions ipi{ x ) = ( n d\'4 > i) an d the hamiltonian 
itself turns into a differential operator. This is due to the fact that two-body interactions conserving 
the 0(5) invariant connect states only with Ana = 0, ±2, which is an infinitesimal change of x. For the 
simplified hamiltonian (56) with x — this procedure yields the following result: 



H osc (r]) 



4 

' N 2 



(i-^)£(i-s 




X — 



(5t?-4) 2 




(97) 



The above hamiltonian is hermitian and can be regarded as a possible quantization of a shifted oscillator 
with a specific x-dependent mass parameter. 

The oscillator potential from Eq. (97) is shown in Fig. 18. The energy Eq and position x of 



the oscillator minimum depend smoothly on the control parameter, as well as the width A l ' 2 . The 
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li^osc (*!>*) 




Figure 18: The oscillator potential from Eq. (97) scaled by (1 — rj) . The physical domain is limited by the 



upper energy E c = where the potential reaches x = — 1. 



departure from the 0(6) limit (77 = 0) is therefore expressed just as a smooth shift and rescaling of 
wave functions and a gradual change of the corresponding energies. However, the solutions loose sense 
if the energy exceeds a certain limit E c , where the semiclassical wave function overflows the physical 
domain x G [— 1,+1]. Fig. 18 shows that E c = (independent of 77), i.e. it coincides with the top 
of the central maximum of the sombrero potential in the present parametrization. At this point, the 
solution becomes singular since the mass diverges and the analytic continuation from the 0(6) limit 
cannot be maintained any more. Therefore, the energy E c = is the upper limit where for each given 
level with r = the concept of quasi-0(6) dynamical symmetry can be applied. Above this energy, the 
level converts to the judicature of another type of analytic continuation that can be named quasi-U(5) 
(cf. Subsec. 4.4). As seen in Fig. 18, the quasi-0(6) energy region (below E c ) shrinks with increasing 



77 and disappears at r) c , where the ground state reaches the critical energy. 

In fact, the same conclusions can be obtained also directly from the classical hamiltonian (96). For 
T rot = and the relative seniority -^ — > (hence p 7 = 0), the sombrero potential creates a phase 
space separatrix at total absolute energy E c = 0. The action along a periodic orbit at energy E can be 
expressed as 



S(v,E) = 2 



0i(v,E) 



Mv,e) 



2M{ m /3)[E-V{ m /3)]dp : 



(91 



where M(r};(3) is a position-dependent "mass" that can be extracted from Eq. (96), while (3q and (5\ 
represent classical turning points of the radial motion. It can be shown that the action as a function 
of energy has a singular tangent at E = E c , when the separatrix is crossed. At this critical energy, the 
inner turning point /3o(jj,E), which for E < E c is a decreasing function of energy, reaches the value 
/3 (77, E) = and gets fixed [100]. 

As discussed by Caprio, Cejnar, and Iachello [117], the singularity at E = has two crucial conse- 
quences: (a) it implies an infinite local growth of the -^ = level density, and (b) creates a singularity 
in the dependence of individual quantum energies on 77. Both these conclusions can be obtained from 
the Bohr-Sommerfeld quantization condition 



S(r), Ei) = 2%h (i + |) = const 



(99) 
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with -J^S equal to the 



dE" 



where K oc N^ 1 . The gap between adjacent levels is given by A ex j^(-^S) 
classical orbit period t p . For E — > E c the period grows asymptotically (it takes an infinite time to get 
over the central maximum of the potential), thus the level density has an anomaly. The variation of 
the ith level energy with 77 can be associated with a contour curve of the function S(j], E) satisfying 
Eq. (99). Exploiting the condition dS = jr-dr/ + M,dE = and evaluating both partial derivatives of 



the action in Eq. (98), one obtains 



d rp 

dr) % 



d S( V ,E) 



dr) 



^S( V ,E) 



E=Ei{rj) 



di) 



«dfo) 



0- 



Pi 



+ W - P 



(100) 



where (/) = j- J Q P f(p, q) dt stands for the time average over the specific periodic orbit. Note that 



Eq. (100) is a semi classical analog of the first Eq. (93). As the level energy crosses the separatrix value 
E c = at some value of the control parameter, r\ = r/ c , the average on the right hand side of Eq. (100) 



drops to zero since the motion gets temporarily frozen at the top of the sombrero potential. A more 
detailed analysis [117] shows that this is a logarithmic type of singularity, with the second derivative 
j-j-Ei approaching to =Foo as 77 converges to r/ c from the left and right, respectively. We therefore 
encounter a continuous quantum phase transition (with no Ehrenfest classification) of excited level i. 

It should be stressed that anomalous properties of wave functions at the classical phase space 
separatrix for a double well potential were studied already in 1987 by Cary, Rusu, and Skodje [245- 
247]. The relation of quantum phase transitions to classical instability points was anticipated by Heiss 
and Miiller [244] and a singular evolution of the quantum spectrum in the above sense has been first 
investigated within the Lipkin model by Leyvraz, Heiss, Scholtz, and Geyer [137,138]. The classical 
limit of the Lipkin model represents a ID projection of the IBM classical hamiltonian for L = and 
jr = 0. If the above analysis within the IBM is extended to states with -^ 7^ (hence p 7 7^ 0), one finds 



out that the centrifugal barrier in Eq. (96) removes the phase space separatrix and the nonanalytic 
behavior of Ei is washed out. Therefore, the ESQPT behavior in the N — > 00 limit applies only to 



states with 



N 



[100]. 



8 Concluding remarks 

In this review, we attempted to summarize the present state of knowledge on quantum phase transitions 
in the interacting boson model. This subject is important from both experimental and theoretical 
viewpoints. It has direct applications in nuclear structure physics, where shape phases and critical 
regions in the chart of nuclides have been subject to intense research in the last decade. Closely related 
models, with analogous types of quantum phase transitions, are used also in molecular physics. 

With respect to recent reviews [23, 24] of the experimental aspects of the structural evolution in 
nuclei, here a stronger accent was put to theoretical problems. Our aim was to present the family of 
interacting boson models as a useful framework for probing the origins and fundamental properties of 
quantum phase transitions in many-body systems, which at the same time holds a close relation to 
experimental data. Rather than going into deep details within a restricted area, we decided to show 
a vast variety of diverse approaches and related problems. The reader interested in a specific issue is 
encouraged to follow the relevant references. 

Basic properties of the interacting boson model connected with phase transitions are contained in 
its simplest version, the IBM-1, which has been mostly discussed in the QPT context so far. One of the 
main advantages of this model is a simultaneous occurrence of first- and second-order quantum phase 
transitions. We showed the roots of the ground-state critical behavior in the mean-field dynamics and 
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outlined the related approaches going beyond the mean field and to the direction of statistical physics. 
These concepts are essential for the description and fundamental understanding of the ground-state 
evolution. 

Various generalizations of symmetry were presented as important guides for the QPT physics of 
many-body systems. The symmetries emerge (somewhat surprisingly) in transitions between different 
dynamical regimes of the system, giving rise to specific signatures of quantum critical behaviors in the 
spectra of low-lying states. The symmetry related aspects of QPT's in many-body systems constitute 
an important subject of ongoing research. 

Extensions of the IBM-1 to various sides, also reported in this review, provide enriched phase 
structures and modified realizations of the QPT phenomena. Phase transitions in these generalized 
models represent an interesting subject on its own, but they can also be seen as concrete cases of a 
general branch of problems concerning the QPTs influenced by a coupling of the system to additional 
degrees of freedom. In the IBM extensions, the main effort has been spent so far to identify the basic 
form of the relevant phase diagrams, while various sophisticated approaches developed in the simpler 
cases still wait for application. 

The last part of the review was devoted to excited-state quantum phase transitions, which represent a 
promising topic for future analyses. In these transitions, global features of the whole quantum spectrum 
are subject to nonanalytic changes of the phase transitional type. New hints for traditional physics of 
phase transitions may follow from the study of these issues. An open problem remains the extension of 
the concept of excited-state QPTs to first-order transitions. The interacting boson model, which itself 
seems to be at the first-order phase transition between simplicity and complexity (with both aspects 
coexisting within one framework), may serve as a firm ground for such studies. 
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